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operations  are  presented.  One  approach  realizes  the  halftone 
method  of  nonlinear  optical  processing  in  real  time  by  replacing 
the  conventional  photographic  recording  medium  with  a real-time 
image  transducer.  In  the  second  approach  halftoning  is 
eliminated  and  the  real-time  device  is  used  directly.  In  this 
case  the  nonlinearity  is  obtained  by  the  inherent  nonlinear 
characteristic  of  the  real-time  device.  A third  method  uses  the 
ability  of  certain  real-time  devices  to  perform  an  intensity- 
to-spatial  frequency  conversion 

Following  a review  of  halftoning  and  real-time  devices,  a 
general  analysis  of  the  halftone  process  considering  the  nonideal 
characteristic  of  the  recording  medium  is  presented.  From  this 
analysis  one  can  predict  the  amount  of  degradation  of  the  output 
due  to  different  parts  of  the  recording  medium  characteristic 
curve  for  any  nonlinear  transformation.  Specific  results  for 
logarithmic  and  level  slice  processes  are  obtained. 

The  problem  of  designing  the  halftone  screen  cell  shape  is 
considered  for  a piecewise  linear  recording  medium  by  solving  a 
nonlinear  integral  equation.  It  is  shown  that  the  solution  can 
be  achieved  for  certain  monotonic  functions  including  exponential 
and  power  transformations.  To  obtain  the  solution  in  general,  an 
approximate  method  which  considers  a discrete  halftone  screen 
density  profile  is  described.  This  gives  the  halftone  screen 
density  profile  for  any  form  of  recording  medium  characteristic 
curve  and  any  type  of  nonlinearity  by  minimizing  in  mean-square 
sense  the  difference  between  desired  and  degraded  outputs. 

The  results  of  computer  simulation  for  logarithmic  and  level 
slice  functions  are  shown.  Experimental  results  are  obtained 
for  an  optimized  halftone  screen  which  matches  the  characteristic 
of  a liquid  crystal  light  valve  (LCLV) . An  overall  logarithmic 
transfer  function  is  produced  in  real  time. 

The  procedure  for  obtaining  nonlinearities  without  halftone 
preprocessing  is  explained.  Experimental  results  using  a 
special  birefringent  liquid  crystal  device  which  produces  a 3-bit 
parallel  A/D  converter  in  real  time  are  shown. 

The  variable  grating  mode  (VGM)  of  liquid  crystal  light 
valve  systems  and  its  application  in  real-time  nonlinear  optical 
processing  is  discussed.  A detailed  analysis  of  the  variable 
phase  grating  is  made  and  from  this  the  processing  limitations 
of  the  VGM  are  determined.  Experimental  results  showing  real- 
time level-slicing  with  a VGM  liquid  crystal  device  are  given. 

Possible  extensions  and  limitations  of  this  work  are 
discussed. 
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ABSTRACT 


Existing  methods  for  nonlinear  optical  processing 
involve  an  intermediate  photographic  step  which  prevents 
real-time  operation.  Although  optical  information 
processors  have  a large  space-bandwidth  product,  the 
problem  of  real-time  input  has  been  a major  limitation  in 
taking  advantage  of  this  processing  ability.  In  this 
research,  the  realization  of  nonlinear  optical  processing 
in  real  time  is  investigated.  Several  techniques  for 
performing  such  operations  are  presented.  One  approach 
realizes  the  halftone  method  of  nonlinear  optical 
processing  in  real  time  by  replacing  the  conventional 
photographic  recording  medium  with  a real-time  image 
transducer.  In  the  second  approach  halftoning  is 
eliminated  and  the  real-time  device  is  used  directly.  In 
this  case  the  nonlinearity  is  obtained  by  the  inherent 
nonlinear  characteristic  of  the  real-time  device.  A third 
method  uses  the  ability  of  certain  real-time  devices  to 
perform  an  intensity-to-spatial  frequency  conversion. 

Following  a review  of  halftoning  and  real-time 
devices,  a general  analysis  of  the  halftone  process 

considering  the  nonideal  characteristic  of  the  recording 
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| medium  is  presented.  From  this  analysis  one  can  predict 
| the  amount  of  degradation  of  the  output  due  to  different 
i parts  of  the  recording  medium  characteristic  curve  for  any 
; nonlinear  transformation.  Specific  results  for  logarithmicj 
I and  level  slice  processes  are  obtained. 
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The  problem  of  designing  the  halftone  screen  cell 

1 

shape  is  considered  for  a piecewise  linear  recording  medium] 
by  solving  a nonlinear  integral  equation.  It  is  shown  that 
the  solution  can  be  achieved  for  certain  monotonic 
1 functions  including  exponential  and  power  transformations. 

! To  obtain  the  solution  in  general,  an  approximate  method! 

| i 

which  considers  a discrete  halftone  screen  density  profile 

I 

] is  described.  This  gives  the  halftone  screen  density 

profile  for  any  form  of  recording  medium  characteristic 

I _ ! 

| curve  and  any  type  of  nonlinearity  by  minimizing  in 

| I 

mean-square  sense  the  difference  between  desired  and 

degraded  outputs.  The  results  of  computer  simulation  for] 

logarithmic  and  level  slice  functions  are  shown. j 

Experimental  results  are  obtained  for  an  optimized  halftone 

i 

i 

screen  which  matches  the  characteristic  of  a liquid  crystal 
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! light  valve  (LCLV) . An  overall  logarithmic  transfer 
] function  is  produced  in  real  time. 

The  procedure  for  obtaining  nonlinearities  without 
halftone  preprocessing  is  explained.  Experimental  results 
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using  a special  birefringent  liquid  crystal  device  which 
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produces  a 3-bit  parallel  A/D  converter  in  real  time  are 
shown. 

The  variable  grating  mode  (VGM)  of  liquid  crystal 
light  valve  systems  and  its  application  in  real-time 
nonlinear  optical  processing  is  discussed.  A detailed 
analysis  of  the  variable  phase  grating  is  made  and  from 

I 

this  the  processing  limitations  of  the  VGM  are  determined. 
Experimental  results  showing  real-time  level-slicing  with  a 
VGM  liquid  crystal  device  are  given. 
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Possible  extensions  and  limitations  of  this  work  are1 
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11-5]  obtained  a logarithmic  transformation  using  this 
method.  Later  Sawchuk  and  Dashiell  [1-6]  performed  a level 
slice  transformation  with  the  halftone  method.  Lohmann  and 
Strand  [1-7]  used  this  method  to  obtain  analog-to-digital 
conversion.  Pseudocolor  with  the  halftone  method  was  shown! 
by  Goodman  and  Liu  [1-8]. 

None  of  the  above  nonlinear  operations  has  been; 

obtained  in  real  time.  This  means  that  the  nonlinear { 

transformation  of  an  input  scene  cannot  be  done  without 
significant  time  delay.  This  is  due  to  the  fact  that  the 

j 

input  and  output  of  these  processors  often  rely  on 

i photographic  techniques  which  are  not  a real-time 

realization.  Recent  developments  have  simplified  the 
output  problem:  television  and  solid  state  devices  are 
available  to  efficiently  make  use  of  the  two-dimensional 

i 

. processed  output.  In  many  situations,  the  human  eye  or 
i observer  is  the  end  user  of  the  information,  so  that 


optical  systems  with  their  inherent  two-dimensional  nature; 


are  ideally  suited 

to 

process  pictorial 

information 

intended  eventually 

for 

the 

human  observer. 

The  major 

difficulty  lies  with 

real- 

-time 

input  to  these 

processors. 

Flexible  real-time 

optical 

input  modulators 

which  can 

convert  image  information 

into 

a form  for  input  to  the 

processing  system  are  needed.  Significant  research  over 
the  last  several  years  has  made  progress  on  certain  aspects 


of  this  problem.  For  ordinary  linear  optical  processing  a 
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great  deal  of  progress  has  been  made  in  developing  fast,! 
sensitive  and  easily  usable  real-time  controllable 
replacements  for  photographic  film.  There  are  many 
different  materials  and  systems,  each  with  its  own 
particular  characteristics  and  tradeoffs  [1-9  to  1-13]  - 
One  common  characteristic  of  nearly  all  current  devices  is 
that  they  are  intended  to  operate  linearly  with  moderate 
contrast  (photographic  gamma  of  3)  over  a, broad  dynamic 
range  (densities  of  0 to  2 or  3).. 


Our  main  goal  in  this  research  is  to  study  the 
feasibility  of  doing  nonlinear  optical  processing  using  a 
real-time  image  converter  and  present  methods  which 
successfully  perform  such  operations.  We  will  take  two 
different  approaches  to  this  problem:  1)  realization  of  the 

i 

halftone  process  in  real  time;  and  2)  obtaining 
nonlinearities  without  halftone  preprocessing.  In  the 
following  we  first  introduce  these  two  methods  and  then 
I present  the  organization  of  this  dissertation  and  its 

i 

research  contributions. 

1.2  Real-Time  Nonlinear  Processing  with  Halftone  Screens 
i ’ ’ 

In  halftone  nonlinear  processing  the  continuous  tone 
input  picture  is  transformed  into  a binary  picture  by 
contact  printing  the  continuous  input  data  through  a 
halftone  screen  onto  a high-contrast  recording  medium.  The 
product  of  the  input  and  halftone  screen  transmittances  is 


clipped  in  the  process,  giving  an  array  of  binary  dots 
whose  size  is  a function  of  clip  level,  the  input 
transmittance,  and  the  halftone  transmittance  profile.  The 
periodic  nature  of  the  halftone  screen  causes  each 
subregion  of  the  binary  image  corresponding  to  a constant 
input  intensity  to  become  quasi-per iodic . When  placed  in 
the  usual  coherent  optical  filtering  system,  multiple 
diffraction  orders  appear  in  the  Fourier  transform  plane 
because  of  the  sampled  input.  The  procedure  for  producing 
nonlinearities  involves  the  use  of  one  diffraction  order 
combined  with  specially  made  halftone  screens.  A filter  is 
placed  in  the  Fourier  plane  that  transmits  the  light  around 
one  diffraction  order  and  blocks  everything  else.  This  in 
effect  demodulates  the  image  [1-3].  After  the  filtered 
diffraction  order  in  the  Fourier  plane  is  inverse 
transformed,  the  continuous  nonlinearly  transformed  output 
appears. 

To  realize  this  process  in  real  time  the  photographic 
hard-clipping  step  should  be  replaced  with  some  real-time 
method.  This  means  that  a real-time  image  converter  is 
needed  instead  of  the  photographic  film.  Among  presently 
developed  devices,  the  liquid  crystal  light  valve  (LCLV) 
appears  to  be  one  of  the  most  promising  for  use  in  optical 
processing  systems  [1-12,1-13],  It  offers  the  advantages 
of  1)  relatively  good  modulation  transfer  function  (to  50 
or  60  cycles/mm.);  2)  fast,  reusable,  simple  operation;  3) 
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room  temperature  operation;  and  4)  physical  compactness 
with  minimal  controlling  electronics  required.  This  is  the 
type  of  device  that  has  been  chosen  for  use  in  real-time 
halftone  processing  experiments. 

Like  most  other  real-time  devices,  the  LCLV  has 
generally  been  optimized  for  a large  linear  dynamic  range, 
in  contrast  to  the  hard-clipping  characteristic  required 
for  halftone  processing.  Given  the  choice  of  a linear  LCLV 
for  real-time  halftone  input,  methods  of  predicting  the 
degradation  and  compensating  for  the  nonideal 
characteristic  have  been  studied.  Dashiell  and  Sawchu. 
have  investigated  some  implications  of  non-ideal  recording 
materials  [1-14],  They  found  methods  to  compensate  for 
finite  gamma  and  finite  saturation  of  the  photographic  film 
for  monotonic  nonlinearities. 

We  will  consider -a  recording  medium  with  a general 

shape  of  characteristic  curve  and  try  to  compensate  for  its 

% 

nonideal  characteristic  for  general  nonlinear  functions. 

1 . 3 Real-Time'  Nonlinear  Processing  Without  Halftone 
Screens 

We  present  two  different  methods  for  obtaining 
noniinearities  without  halftone  screens;  a)  direct;  b) 
variable  grating  mode. 


1.3.1  Direct  Nonlinear  Processing 


To  achieve  nonlinear  functions  without  halftone! 
i screens  the  inherent  nonlinear  characteristic  of  a 

i 

recording  medium  or  real-time  transducer  is  directly  used. 
This  invovles  the  proper  biasing  and  selection  of  operating 

points  on  a nonlinear  curve  to  directly  achieve  a 

point-by-point  intensity  transformation.  The  processing 

i 

may  take  place  with  incoherent  illumination,  avoiding  the 

! 

problems  of  speckle,  phase  noise,  and  the  possible 

necessity  for  a laser  source.  It  also  requires  a lower 
space-bandwidth  product  on  the  real-time  input  device 
; because  the  halftone  screen  is  eliminated.  One  drawback  is 
that  the  flexibility  of  this  technique  is  considerably  more 
limited  than  halftoning.  Tai,  Cheng  and  Yu  [1-15]  have 

I 

obtained  a logarithmic  nonlinearity  directly  (not  in  real 
time)  using  photographic  film  with  special  processing.  . 


1.3.2  Variable  Grating  Mode  Nonlinearities 


This  method  involves  the  conversion  of  different  input | 
intensities  to  a local  phase-modulated  grating  whose! 
spatial  frequency  is  a function  of  the  brightness.  When 
the  variable  grating  is  placed  in  the  front  focal  plane  of 
a coherent  Fourier  transform  processor,  the  difference  in 
local  spatial  frequency  should  cause  different  input  levels 
to  be  effectively  placed  at  different  points  in  the 
transform  plane.  By  selective  filtering  and  recombination 

j 

of  the  transform  components,  various  nonlinearities  should 


be  possible.  This  method  relies  on  the  behavior  of  certain 
| liquid  crystal  real-time  devices  which  have  shown  the 

I 

I variable  grating  mode  in  experiments  [1-16]. 

j 1.4  Organization  of  the  Dissertation  and  Its  Research 

Contributions 
■ — ■ — ■"  * ■'■■■■ 

I 

I ; 

In  the  second  chapter  we  review  halftone  nonlinear 

I , 

| processing  and  characteristics  of  a liquid  crystal  light 

! valve  (LCLV)  as  a real-time  image  converter.  We  then 

| 

explain  how  these  two  techniques  can  be  combined  to  produce 
, nonlinearities  in  real  time.  Chapter  3 presents  a 

formulation  of  the  halftone  process  with  a nonbinary 
recording  medium.  This  formulation  relates  the  input  and 
output  intensities  through  the  halftone  screen  parameters 
and  the  recording  medium  characteristic  curve.  The  results 
for  the  binary  recording  medium  are  obtained  as  a special 

- 

case  of  the  general  recording  medium. 

In  Chapter  4 we  use  the  formulation  of  chapter  three j 
to  predict  degradations  due  to  the  nonbinary  characteristic 
of  the  recording  medium.  The  result  of  computer 
simulations  for  different  parameters  of  the  recording 
medium  characteristic  are  shown. 

In  Chapter  5 we  try  to  compensate  for  the  degradations 

| 

by  analytical  methods.  It  is  shown  that  for  some 

I 

j nonlinear ities  we  can  obtain  the  halftone  screen  parameters 


without  any  precompensation  for  the  nonbinary  recording 
medium. 

Chapter  b gives  a new  analysis  of  the  halftone  process 
in  which  the  halftone  screen  density  function  is  considered 
as  a discrete  function.  With  such  analysis  we  describe  new 
optimized  methods  of  halftone  screen  design.  Results  of 
computer  simulation  for  different  parameters  of  the 
recording  medium  characteristic  are  shown. 

In  Chapter  7 we  discuss  methods  of  making  halftone 
screens  and  present  experimental  results  for  a real-time 
logarithmic  process  using  an  optimized  halftone  screen  with 
a LCLV. 

Chapter  b explains  how  the  nonlinear  characteristic  of 
a special  birefringent  LCLV  can  be  used  to  generate  an 
analog-to-dig ital  conversion  nonlinearity  directly  without 
any  halftone  preprocessing.  The  result  of  an  experiment 
producing  a 3-bit  parallel  A/D  conversion  in  real  time  is 
shown. 

In  Chapter  9 the  variable  grating  mode  VGM  of  the 
liquid  crystal  and  its  application  in  real-time  nonlinear 
processing  is  discussed.  In  a detailed  analysis  we 
determine  the  processing  limitations  of  the  VGM  effect. 
Preliminary  experimental  results  using  a VGM  liquid  crystal 
device  for  real-time  level  slicing  is  presented. 
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Chapter  10  discusses  conclusions  and  possible' 
extensions  of  this  work. 

The  specific  research  contributions  of  this  work  are 
now  described.  The  analysis  of  the  halftone  nonlinear 
process  has  been  extended  to  predict  the  nonlinear 
input-output  curve  with  any  recording  medium  characteristic 
curve  shape  and  any  halftone  screen  cell  shape.  New 
analytical  design  formulas  for  halftone  cell  shape  with  a 
piecewise  linear  recording  medium  and  certain 
nonlinearities  have  been  derived.  New  design  procedures 
for  halftone  screen  cell  shape  are  presented  which  produce 
screens  that  compensate  for  the  nonideal  characteristic  of 
the  recording  medium  with  any  nonlinearity.  New  methods 
have  been  developed  for  obtaining  real-time  nonlinear 
functions  without  halftone  preprocessing  using  the 
birefringent  liquid  crystal  device  and  the  variable  grating 
mode  (VGM)  liquid  crystal  device.  Experimental 
demonstrations  of  logarithm,  3-bit  parallel  A/D  conversion, 
and  level  slice  functions  in  real  time  have  been  made. 


CHAPTER  2 


REVIEW  OF  HALFTONE  NONLINEAR  PROCESSING  AND 
THE  LIQUID  CRYSTAL  LIGHT  VALVE 

In  this  chapter  we  review  nonlinear  halftone 
processing  and  the  liquid  crystal  light  valve  (LCLV)  as  a 

| 

I real-time  image  converter.  We  then  discuss  how  they  can  be 
' combined  to  produce  parallel  nonlinear  transformations  in 
real  time. 

2. 1 Nonl inear  Optical  Processing  with  Halftone  Screens 

I 

Halftone  nonlinear  optical  processing  methods  rely  on 
!a  pulse  -width  modulation  intermediate  step,  in  which  the 
input  data  is  coded,  and  a filtering  part  that  recovers  a 
i continuous-tone  picture  from  the  coded  one  [2-1],  The i 
pulse-width  modulation  step,  often  called  the  halftone 
process,  transforms  a continuous-tone  picture  into  a 
l picture  with  only  two  levels  of  density.  The  resulting 

I 

I picture  is  called  a halftoned  picture.  It  consists  of 

I 

groups  of  dots  or  bars  depending  on  whether  the  halftone 
screen  used  in  the  process  is  one  or  two-dimensional. 
Corresponding  to  each  density  level  in  any  part  of  the 
original  picture  a group  of  equal-area  dots  or  a series  of 
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equal-width  bars  is  present  in  the  halftoned  picture.  The 
area  of  the  dots  or  the  width  of  the  bars  is  a function  of 
the  corresponding  density  level  in  the  original  picture  and 
the  density  profile  of  the  halftone  screen.  Thus,  the 
information  about  the  original  picture  density  levels  is 
encoded  into  areas  of  the  dot  groups  or  the  widths  of  bars 
present  in  the  halftoned  picture.  The  method  is  shown  in 
Fig.  2.1.  The  input  illumination  passes  through  the 
original  picture  transparency  and  falls  on  a periodic  array 
of  partially  absorbing  dots  or  lines  (halftone  screen). 
The  light  transmitted  by  the  halftone  screen  falls  upon  a 
hard-clipping  photograpic  film  which  acts  as  a thresholding 
device  and  records  a binary  image.  The  operation  for  a one 
dimensional  halftone  screen  is  shown  in  Fig.  2.2.  Note 
that  the  densities  of  the  input  picture  transparency  . and 
the  halftone  screen  add  together  to  effectively  partially 
absorb  this  illumination  at  any  point.  Depending  on 
whether  this  sum  of  densities  is  less  than  or  greater  than( 
the  threshold  density  of  the  hard-clipping  film,  it  will  or 
will  not  be  activated.  This  causes  the  amplitude 
transmittance  for  the  encoded  picture  transparency  to  look 
like  Fig.  2.2(e).  From  the  amplitude  transmittance  of  the 
halftoned  picture  we  see  that  the  effect  of  the  halftone 
process  on  the  input  picture  is  equivalent  to  sampling.  If 
this  sampled  picture  is  used  as  an  input  to  a coherent 
optical  processing  system  as  shown  in  Fig.  2.3,  different 
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Figure  2.1  Production  of  half toned  picture 


diffraction  orders  appear  in  the  filter  plane.  To  be  able 

| 

to  successfully  filter  out  any  single  order,  the  spatial 
frequency  of  the  halftone  screen  should  be  much  greater 
than  the  maximum  spatial  frequency  present  in  the  input 
picture.  Assuming  that  such  a requirement  holds,  we  can 
put  a simple  spatial  filter  in  the  Four ier.  plane  and  select 
one  of  the  orders.  The  selected  order  is  then  transformed 
back  to  give  a desampled  image  in  the  output  plane  of  the 
processor.  The  intensity  at  any  point  of  this  output 
picture  is  then  a nonlinear  function  of  the  corresponding 
point  in  the  input.  The  form  of  this  nonlinearity  is 
determined  by  the  shape  of  the  halftone  screen  density 
profile  and  the  diffraction  order  selected  in  the  spatial 
filtering  system.  Details  of  the  mathematical  relationship 
between  the  output  and  input  intensities  will  be  worked  out 
in  the  following  chapters. 

2 . 2 Liquid  Crystal  Light  Valve 

The  term  "liquid  crystal"  is  applied  to  substances 
whose  rheological  behavior  is  similar  to  that  of  fluids  but 
whose  optical  behavior  is  similar  to  the  crystalline  state 
over  a given  temperature  range  [2-2].  One  class  of  liquid 
crystals  called  nematics  is  the  most  common  type  used  in 
devices  whose  properties  interest  us.  The  nematic  liquid 
crystal  phase  is  formed  by  many  organic  compounds.  The 
most  significant  common  feature  of  all  these  compounds  is 
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that  they  have  long  rodlike  molecules  [2-3].  When  an' 

] 

electrical  field  is  applied  to  a liquid  crystal,  the 

! * * j 

material  is  reoriented  in  space  by  both  the  field  and  ionic 
conduction  effects.  If  the  conductivity  of  the  liquid! 
crystal  is  sufficient,  conduction  effects  will  predominate, 
whereas  in  the  absence  of  significant  conductivity  effects, 
the  nematic  is  simply  reoriented  by  the  field  [2-4]. 

The  field  effect  operation  of  a LC  device  is  based  on 
'the  electrically  controlled  birefringence  characteristic  of 
the  nematic  liquid  crystal.  The  LC  layer  is  normally 

placed  between  a transmitting  electrode  and  a reflecting 
mirror.  An  incident  polarized  beam  of  light  on  the 
transmitting  electrode  will  be  reflected  back  from  the 
mirror  while  traveling  through  t^he  liquid  crystal  as  shown 
in  Fig.  2.4. 

Depending  on  the  voltage  across  the  liquid  crystal, 

' 

the  reflected  light  will  experience  some  phase  retardation! 
due  to  reorientation  of  the  LC  molecules  by  the  field. 
This  phase  retardation  can  be  detected  as  an  intensity 
variation  if  we  pass  the  reflected  light  through  a second 
polarizer  that  is  either  parallel  or  perpendicular  to  the 
one  in  the  incident  beam.  With  such  arrangement  we  have  a 
light  valve  in  which  the  electric  field  variation  across 
the  LC  controls  the  intensity  variation  of  the  output  beam 
(reflected  light).  Now,  if  we  introduce  a photoconductor 
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Figure  2.4 


Diagram  of  a tunable  birefringence  field 
effect  cell  (P=Polarizer , A=Analyzer, 
TE=Transparent  electrode,  M=Mirror, 
LC=Liquid  crystal 
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layer  in  the  system  we  can  transform  the  intensity 

variation  of  a second  light  source  (input  light)  to  voltage 

variation  across  the  photoconductor  that  in  turn  can  be 

transferred  to  voltage  variation  across  LC.  Hence,  the  end 

result  is  a light  controlled  light  valve.  This  is  shown  in 

Fig.  2.5  with  a light  blocking  layer  added  to  separate  the 

o 

input  and  output  beams.  If  the  LC  is  twisted  by  45  we 
could  get  a dark  off-state  and  the  resulting  device  would 
be  the  so-called  standard  linear  LCLV  [2-5].  This  device 
can  be  made  to  operate  with  a speed  of  about  30 
frames/second  or  TV  frame  rate. 

2 . 3 Real-Time  Nonlinear  Optical  Processing  with  Halftone 
Screens 

To  realize  the  halftone  process  in  real  time  both  the 
halftoning  step  and  output  readout  should  be  performed 
rapidly.  As  described  in  chapters  1 and  2 the  output  is 
easy  to  obtain  in  real  time.  To  achieve  the  halftoning  in 
real  time  the  photographic  hard-clipping  medium  should  be 
replaced  with  a real-time  image  converter.  The  ideal 
real-time  device  should  have  two  properties:  1)  the  ability 
to  convert  an  incoherent  image  to  a coherent  image;  2) 
infinite  gamma.  Everywhere  in  this  work  we  define  the 
parameter  gamma  as  the  slope  of  the  linear  portion  of  the 
output  amplitude  transmittance  vs.  log  of  input  intensity 
of  the  real-time  image  transducer.  Our  definition  of  gamma 
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Figure  2.5 


Cross-sectional  schematic  of  LCLV 


should  not  be  confused  with  the  photographic  definition  of 
gamma  as  the  slope  of  the  linear  part  of  density  vs.  log  of 
the  input  intensity.  The  first  property  is  common  to  many 
image  transducers  and  the  liquid  crystal  light  valve  (LCLV) 
described  shares  this  property.  On  the  other  hand,  almost 
all  real-time  image  converters  are  optimized  for  linear 
processing  and  do  not  possess  the  second  property.  This 
can  be  seen  from  Fig.  2.b  where  a typical  characteristic  of 
a standard  linear  LCLV  is  shown. 

In  the  following  chapters  we  investigate  this  problem 
by  developing  a new  formulation  of  the  halftone  process  and 
presenting  methods  for  the  compensating  for  nonideal 
characteristic  of  real-time  image  converters. 


I 
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CHAPTER  3 

HALFTONE  PROCESSES  WITH  GENERAL  RECORDING  MEDIUMS 

Halftone  nonlinear  processes,  described  in  chapter 
two,  have  been  formulated  with  the  assumption  that  a binary 
recording  medium  is  used  in  the  halftoning  step  [3-1]. 
With  this  assumption,  once  the  output  intensity  is 
expressed  as  a function  of  input  intensity  and  the  halftone 
screen  density  profile  (analysis) , we  can  easily  invert  the 
problem  and  get  the  halftone  screen  density  profile  given 
the  relationship  between  the  output  and  input  intensities 
(synthesis) . Unfortunately,  almost  all  recording  media 
deviate  from  the  binary  assumption.  This  deviation  which 
is  quite  small  for  high  contrast  photographic  films,  is 
quite  noticeable  for  any  real-time  image  converter  at  the 
present  time  [ 3-2 , 3-3 ] . Consequently,  to  appreciate  the 
halftone  processes  fully,  we  should  remove  the  assumption 
of  a binary  recording  medium  from  the  very  beginning  of  our 
mathematical  formulation  of  these  processes.  Dashiell  and 
Sawchuk  have  considered  the  effects  of  finite  gamma  and 
saturation  density  of  the  recording  medium  on  halftone 
process  by  modifying  the  formulas  for  the  ideal  recording 
medium  [3-4]. 
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This  formulation  does  not  predict  the  effects  of  the 
curved  portions  of  the  recording  medium  characteristic 
curve  on  the  overall  nonlinear  transformation.  Moreover 
the  formulation  is  restricted  to  monotonic  halftone  cell 
shapes.  In  this  chapter  we  present  a formulation  of  the 
halftone  process  which  considers  a recording  medium  with  a 
characteristic  curve  of  general  shape  and  predict  the  final 
degradation  of  the  output  due  to  such  characteristic  curve 
for  any  halftone  screen  cell  shape.  This  formulation  is 
used  for  the  case  of  a binary  recording  medium  and  its 
result  is  compared  with  previous  derivations. 

3.1!  Mathematical  Formulation 

I 

| For  a general  recording  medium  in  the  halftone  process 
we  should  replace  the  binary  characteristic  curve  of 
Fig.  2-2a  by  a curve  like  Fig.  3.1.  This  causes  the 
amplitude  transmittance  of  the  halftoned  picture  to 
consist  of  pulses  that  are  no  longer  rectangular  as  shown 
in  Fig.  3.2.  The  amplitude,  width  and  shape  of  these 
pulses  depend  on  the  input  picture  density  levels,  halftone 
screen  density  profile,  and  the  shape  of  the  characteristic 
curve  of  the  recording  medium.  Each  group  of  these  pulses 
corresponds  to  a constant  intensity  subregion  in  the  input 
picture.  The  period  L (Fig.  3.2)  of  the  halftone  screen  is 
chosen  to  be  small  in  comparison  with  the  period  of  the 
highest  spatial  frequency  component  in  the  input  picture, 
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L 


so  any  local  region  of  the  amplitude  transmittance  of  the 
halftoned  transparency  is  approximately  a periodic  sequence 
of  pulses  and  can  be  expanded  in  complex  Fourier  series 

+ 00 

th(x)  = ^ B.exp  (-j2irkx/L)  (3.1) 

k=-°° 


where 


*k 


j2irkx 
J T 

t^(x)e  dx. 


(3.2) 


In  the  above  sum  each  term  represents  a grating  order,  and 
when  we  take  the  Fourier  transform  of  the  halftoned  picture 
in  the  coherent  optical  processor,  these  orders  appear  in 
the  Fourier  plane  as  isolated  spectral  islands.  The 
spatial  filter  in  this  plane  selects  a single  order. 
Hence,  the  resulting  intensity  distribution  at  the 
processor  output  is 


1 5k  I 


■'u: 


. 2Tikx 


th (x) e 


dx  | 


(3.3) 


which  relates  the  intensity  at  any  point  of  the  output 
picture  to  the  amplitude  transmittance  of  the  halftoned 
picture  and  the  selected  order.  Now  th(x)  can  be  related 
to  the  input  intensity  as  follows.  Let  the  local  input 
picture  intensity  that  produced  the  above  train  of  pulses 


! 


l 
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If  the  density 


on  the  halftoned  picture  be  denoted  by  I^n 

variation  of  one  period  of  the  halftone  screen  is 

represented  by  f(x),  then  the  intensity  transmitted  by  the 

— f ( x ) 

halftone  screen  is  Iin*10  [3-5],  If  we  let  the 

amplitude  transmittance  versus  log  exposure  curve  of  the 
recording  medium  be  described  by  g(logE),  then  we  can  write 

th(x)  = g{log[Iin*10  = g [logIin~f (x) ] (3.4) 

and  replacing  it  in  Eq.  (3.3),  we  have 

L j2TTkX 

Io(Iin'k)  = f g(logIin-f (x)]e  L dx|2  (3.5a) 

which  relates  the  intensity  at  any  point  of  the  output 
picture  to  the  intensity  of  the  corresponding  point  in  the 
input  picture  through  a nonlinear  integral  relationship. 
When  the  specific  forms  of  g(logE)  and  f(x)  are  substituted 
in  this  relationship  and  the  integral  is  solved,  the 
overall  relation  between  the  output  intensity  and  the  input 
intensity  is  nonlinear.  This  nonlinearity  depends  on 
g(logE),  f(x)  and  the  value  of  the  order  selected. 

The  above  formulation  could  be  performed  in  terms  of 
the.  intensity  transmittance,  t(x),  of  the  halftone  screen 
rather  than  its  density  profile.  In  this  case,  for 
simplicity,  we  also  need  the  transmittance  versus  exposure 
curve  of  the  recording  medium,  T ( E ) , rather  than  g(logE). 
If  we  go  through  similar  arguments  as  above  we  get 


L j 4.  II  A.  A 

IolIi„’k»  - It  I „ T[Iin-T<x)le  L dx|2  . 

3 . 2 Binary  Recording  Medium 

The  characteristic  curve  of  a binary  recording  medium 
is  shown  in  Fig.  3.3.  Ideally  a=0,  and  b=l.  Note  that 
this  form  of  characteristic  curve  is  applicable  to  a 
positive  transparency.  We  could  also  consider  the  more 
familiar  negative  transparency  although  the  basic  results 
remain  the  same.  We  will  choose  the  positive  transparency 
curve  because  they  are  more  similar  to  the  characteristic 
curves  for  real-time  devices.  We  now  simplify  the  general 
relationship  of  Eq.  (3.5a)  using  the  charac cer istic  curve 
of  Fig.  3.3. 


(3.5b) 


3.2.1  Zero  order 

For  the  zero  order  case  k=0 , and  Eq.  (3.5a)  becomes  i 

i 

L 


I din'0*  = {f 
o in  Xj 


g[logIin-f (x) ]dx>  . 


(3.6) 


Considering  Fig.  3.3  we  can  write 


if  logIin-f(x)  < logIr 


then 


f(logE)  = a (3.7) 


and 


if  logIin~f(x)  > logIr  then  g(logE)  * b.  (3.8) 
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Assuming  f(x)  to  be  a monotonical ly  increasing  function  we 
have 


if  x > f-1(log^) 

r 


then  g(logE)  = a 


(3.9) 


1 

if  x < f (log-^) 

r 


then  g(logE)  = b. 


(3.10) 


When  used  in  Eq.  (3.6),  this  gives 


= dr  f 


f'do^)  o 


Vdn'**  - {Z[  J0  bdx 


r 

+ 

. 


i I • 

--1 in. 
f ( log-j — ) 

r 


adx] }2 (3. 11) 


or  after  simplification 


loUin.W  - [a+i^f^dog^i)]2 


(3.12) 


which  is  the  same  as  the  results  obtained  previously  for 
the  case  a=U,  b=l  [3-1]  except  for  the  fact  that  it  is 
obtained  in  a more  straightforward  manner  and  can  easily  be 
generalized . 


3.2.2  Nonzero  order 


For  k=0 , the  above  simplifications  for  the 
characteristic  curve  can  be  used  in  Eq.  (3.5)  to  obtain 


i 1 • 

f_1(log^-)  j 2tt)cx 

I U.  ,k)  = d|  r be  L dx+ 

o in  LJ 0 j 


dirkx 

3 L . , ,2 
ae  dx] | . (3.13) 


f1,  log  ill) 


let 


x x = f 1(log^) 
xr 


then 


2iTkx, 


Io^Iin'k) 


L -1'!2 


2iTkxn 


2Trkx, 


(b-a) 

2 2 
4ir  k 


-^2  . 3- 


(e 


-1)  (e 


"1- 


(b-a)2  „ __"kxl, 

. 2.  i 1— cos  ^ ) 

4tt  k 


-1) 


(3.14) 


(b-a)2  . 2 *kxl 

5TT"  Sln  — 

ir  k 


after  replacing  for  x^  we  have 


I0(Iin-k)  - sin2tT 


in, 


it  k 


(3.15) 


This  also  agrees  with  the  previously  obtained  results 


[3-1]. 

3.3  Halftone  Screen  Density  Functions  for  Some  Useful 
Nonlinearities 


The  halftone  screen  density  function  for  different 
nonlinearities  with  the  general  recording  medium  will  be 
determined  in  later  chapters.  Here,  for  our  reference  we 
will  find  the  halftone  screen  density  functions  for  the 
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[binary  recording  medium  for  several  specified 
nonlinear  1 ties.'. 

I 

3.3.1  Logarithmic  Transformation 

This  transformation  can  be  obtained  in  the  zero  order. 
Hence  from  Eq.  (3.14) 

i 

Io(Iin'0)  = ta+-^^f-1(log-^)]2.  (3.16) 

For  a logarithmic  transformation 

I 

Io(Iin,0)  = Klo9(Iin/Ir)  (3.17) 

where  K is  a constant.  Combining  Eq.  (3.17)  and 
Eq.  (3. 16)  , we  get 

Klog(^)  = [a+-^i  f_1(log^)]2.  (3.18) 

Xr  L r 

li- 
lt we  let  log —iiibe  represented  by  f(x)  and  note  that 

n Ir 

f_i[f(x)]=x  then  Eq.  (3.18)  can  be  written  as 

f(x)  = l[a+^x]2  (3.19) 

which  gives  the  density  profile  of  the  halftone  screen  for 
a logarithmic  transformation. 

3.3.2  Exponential  Transformation 

This  transformation  can  also  be  obtained  in  the  zero 

order.  We  want 
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V'in'01  ' “(6)  Jr 


(3.20) 


where  a and  3 are  constants.  Combining  Eq.  (3.20)  and 
Eq.  (3.16)  we  get 


io\  f . (b-a)  £-l,.  Iin>12 

a (3)  = [a+  ■ -jj  - - f (log-y— )] 


(3.21) 


I”  log(^) 

c,(8)L10 


[a+i^aif-idog^ln2 

r 


(3.22) 


From  this  we  can  write 


,[lof(X)]  = r.a^.Ul2 


[a+33^x]‘. 


(3.23) 


Taking  the  logarithm  of  both  sides  twice  and  simplifying, 
we  obtain 


f(x)  = log{21og  [a+  ^ ^-tx]  -loga)-log  log3.  (3.24) 

Lt 


This  gives  the  density  profile  of  the  halftone  screen  for 
an  exponential  transformation. 


3.3.3  Power  Transformation 


This  is  another  example  of  a transformation  possible 
in  the  zero  order.  We  want 


i 


I 


1 • , 

Io(Iin'0)  = Y(t~)  (3.25) 

r 

where  y and  X are  constants.  Combining  this  with 

Eq.  (3.16)  we  have 

y(_L^)X=  [a+  (b~-a-?  f-1  (log-j^l)  ) 2 (3.26) 

r 1r 

or 

Xiog^  j. 

Y (10)  r = [a+-(-b-I-^-)-f  1 ( log— ?^)  ] 2 (3.27) 

r 

from  which  we  can  write 

Y(10)Af(x)  = [a+(b-a)x]2^  (3.28) 

Taking  the  logarithm  of  both  sides  of  Eq.  (3.28)  and 
simplifying  we  obtain 

f(x)  = |-log  [a+  — ~-aj  x]  logy  . (3.29) 

This  gives  the  density  profile  of  the  halftone  screen  for  a 
power  transformation. 

3.3.4  Level  Slice  Transformation 

In  this  case  the  desired  input  and  output  intensities 

are 
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° for  Iin  •=  ri 

ro  ’ c for  rl  i ri„  < *2 

0 for  I„  < I. 

2 — xn 


(3.30) 


and  this  equation  is  shown  in  Fig.  3.4.  To  obtain  this 
transformation  we  should  use  the  first  order,  i.e.,  k=l  in 
Eq.  (3.15).  We  have 


V1™'11  - si“2lEf'1,109Ta>l  <3-31> 


Combining  Eq.  (3.30)  and  Eq.  (3.31)  we  can  see  that 


f’1(log^)  = 0 for  Iin  < Ix, 


f'dog^)  ^sin"1^  for  Ix  < x.n  < I,, 


-1  I. 

f (log-j^-)  = L for  Iin  >_  1 2. 
These  are  summarized  in  Fig.  3.5  where 


0.32) 


_ _ L . -1  cn 
r " ? Sln  E=i 


(3.33) 


From  Fig.  3.5  we  can  obtain  the  density  profile  o.f  the 
hal.ftone  screen  for  the  level  slice  transformation  as  shown 
in  Fig.  3.6. 


Hlog  lm/lr) 


0 log(I,/Ir)  log  ( I2/Ir) 


Figure  3.5  Inverse  function  of  the  density  profile  of 

the  halftone  screen  which  performs  the  level 
slice  transformation  of  Figure  3.4 


Figure  3.6  Density  profile  of  the  halftone  screen  which 
performs  the  level  slice  transformation 
of  Figure  3.4 


CHAPTER  4 


DEGRADATION  OF  THE  OUTPUT  INTENSITY  IN  THE 
HALFTONE  PROCESS 

I 

The  nonbinary  characteristic  of  the  recording  medium 
in  the  halftone  process  is  the  main  source  of  degradation 
of  the  output  intensity.  This  analysis  theoretically 
predicts  the  departure  of  the  nonlinearity  achieved  from 
the  ideal  so  that  a compensation  for  this  effect  can  be 
included  in  the  design  of  the  halftone  screen.  A typical 

i 

characteristic  curve  of  a nonbinary  recording  medium  is 
shown  in  Fig.  4.1.  Comparing  this  figure  with  the 
characteristic  curve  of  a binary  medium  shown  in  Fig.  3.3, 
we  can  see  that  there  are  three  distinct  regions  in  these 

two  curves  that  make  them  different.  As  shown  in  Fig.  4.1,| 

| 

they  are  usually  referred  to  as  the  toe,  the  linear  part,! 
and  the  shoulder  of  the  characteristic  curve.  The  effect 

l 

of  finite  gamma  and  finite  saturation  of  the  recording 
medium  on  the  halftone  process  has  been  considered  before 
[4-1].  in  the  following  analysis,  we  will  show  how  each  of 
these  three  regions  of  the  characteristic  curve  contributes 
to  the  final  degradation  of  the  output  intensity.  The 
analysis  will  be  done  separately  for  the  zero  and  nonzero 


orders . 


4. 1 Zero  Order 

The  general  formula  relating  the  output  intensity  in 
the  zero  order  to  the  input  intensity  from  Eq.  (3.6)  is 

1 2 

Io(Iin,0)  = J gtlogIin-f(x)]dx}  . (4.1) 

To  find  out  how  different  parts  of  the  nonbinary  recording 

medium  characteristic  curve  affect  I^(I.:„,0),  we  can  use 

o in 

the  halftone  screen  density  profile  that  is  designed  for 
the  binary  recording  medium  and  some  model  characteristic 


curves  which 

are  binary  in  every 

respect  except 

for 

their 

toe,  linear 

part, 

or  shoulder 

as 

shown  in 

Fig . 

4.2, 

Fig.  4.3,  and  Fig. 

4.4.  This 

has 

been  done 

for 

the 

logarithmic  process  as  explained  in  the  following  sections. 

4.1.1  Degradation  Due  to  the  Linear  Toe  of  the  ‘Recording 
Medium  Characteristic  Curve  in  the  Logarithmic  Process 

In  this  case 


VW01  - K 


(4.2) 


where  K is  a constant  and  Ir  is  shown  in  Fig.  3.3.  For  the 
ideal  recording  medium  the  halftone  screen  density  profile 
from  Eq.  (3.19)  is 

f (x)  = i[a+-^^ix]2.  (4.3) 


■ 


L 
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Let 


t = log(Iin/Ir> 

and  (4.4) 

e = log(I2/I1) 

and  now  use  Eq.  (4.2),  Eq.  (4.3)  and  the  characteristic 
curve  of  Fig.  4.3  in  Eq.  (4.1).  For  simplicity  assume  a=0, 
b=l.  After  some  simplification  the  result  is 


2 

/4-^-£  (t+|)3  for  logl.  < logl2 

9e 


i,2aK1/2  r ,.e.  3/2  e . 3/2 , 1/2  , . w.  £.1/2,2 

i — j-e — [ (t+j)  - (t-j)  ]+K  (1-a)  (t-  2)  } 


Io(Iin'0,-< 


for  logI2  < logIin  < R 


(4.5) 


1/2 


, „ 1 / 2 , x 1/2  , a la_  £,  2otKJ-/  ex3/2  a ^,2 

[K  (l-a)(t-y)  +-(t-T) — (t-7)  -y^+aJ 


3eK 


for  ^ < logIin 


The  quantity  a is  defined  in  Fig.  4.3.  These  equations 
have  been  calculated  for  ot  =0.1  and  different  e's.  This 
means  that  the  height  of  the  toe  is  kept  constant  while  its 
slope  changes.  This  way  the  input  intensity  range  over 
which  the  linear  toe  exists  is  extended  to  higher  values  of 
e.  The  result  is  shown  in  Fig.  4.5.  From  this  figure  it 
is  seen  that  increasing  e causes  the  transfer  curve  to  lose 
some  dynamic  range  while  its  shape  is  preserved.  If  we 
assume  a curved  toe  rather  than  a linear  toe,  we  expect  to 
get  some  small  change  in  the  shape  of  the  transfer  curve 
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but  same  loss  of  dynamic  range. 


4.1.2  Degradation  Due  to  the  Linear  Part  of  the  Recording 
Medium  Characteristic  Curve  in  the  Logarithmic  Process 


Using  the  assumptions  of  section  4.1.1  and  the 
characteristic  curve  of  Fig.  4.4  in  Eq.  (4.1)  we  get  after 
some  simplification 


^V-§>3 

9 e 


for  logl-  < logl. 

J in  . 


Io(Ijn'0)  = <t+f)3/2-(t-f)3/2]2  for  logI2<logIin<i  (4.6) 


1 /x.  £\  2K  /4 — £\  3/2 1. 1 i 2 ~ 1 


3^1 1 for  f 1 logIin 


where  we  assume  a=0  and  b=l.  These  equations  have  .been 
calculated  for  different  values  of  e.  The  result  is  shown 
in  Fig.  4.6.  From  this  figure,  we  can  see  that  extending 
the  linear  part  of  the  recording  medium  characteristic 
curve  to  larger  input  intensity  values  causes  the  transfer 
curve  to  lose  dynamic  range  and  also  to  change  shape. 
Although  the  loss  in  dynamic  range  is  undesirable,  the 
transfer  curve  shape  change  is  even  less  desirable. 


4.1.3  Degradation  Due  to  the  Linear  Shoulder  of  the 
Recording  Medium  Characteristic  Curve  in  the  Logarithmic 


Process 


I 

1 
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Using  the  assumptions  of  section  4.1.1  and  the 
characteristic  curve  of  Fig.  4.5  in  Eq.  (4.1)  we  get  after 
some  simplification 

K[2(3~e)  (t+|)3/2+8(t+|)1/2]2  for  logIin<logI2 

R{2  (3~-B)-[  (t+|)  3/2-(t-|)  3/2]+6(t+|)1/2}2 

for  logI2<logIin<l/K  (4.7) 

r (1-6)  #J.  e,  2(1-6)  ^ £>3/2,2 
3e  {t~T  1 

for  i < logI.n  . 

where  we  assume  a=U  and  b=l  as  before.  The  result  of  the 
above  equations  for  8=.9  and  different  values  of  £ are 
shown  in  Fig.  4.7.  The  transfer  curve  change  is  seen  to  be 
similar  to  the  one  for  the  linear  toe. 

Noting  that  the  input  intensity  range  over  which  the 
toe  and  shoulder  of  the  recording  medium  characteristic 
curve  are  present  is  usually  much  smaller  than  the  range 
for  the  linear  part,  we  can  conclude  that  the  linear  part 
is  responsible  for  most  of  the  degradation  in  the  output 
intensity  for  a monotonic  nonlinearity  like  the  logarithmic 
process. 

4. 2 Nonzero  Order 

The  general  formula  relating  the  output  intensity  in 
the  nonzero  order  to  the  input  intensity  from  Eq.  (2.5)  is 
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. 2irkx 


Xo  (Iin,lc)  = IT  L gUogIin-f  (x)]e  L dx|2 


(4.8) 


' 

i where 


is  the  order  selected  and  the  rest  of  the' 


parameters  are  explained  in  section  4.1.1.  To  find  the 

I 

degradation  due  to  each  part  of  the  characteristic  curve  of 
! the  recording  medium,  we  can  go  through  the  same  sort  of 
: procedures  as  in  the  zero-order  case.  The  specific  example 
we  consider  here  is  the  level  slice  function.  It  turns  out 
in  this  case  that  the  calculations  simpl ify  and  thedif ferent 
parts  of  the  recording  medium  characteristic  curve  do  not 
have  to  be  considered  separately.  This  is  shown  in  the 
following  sections. 

4.2.1  Degradation  of  the  Output  Intensity  Due  to  the 
Nonbinary  Characteristic  of  the  Recording  Medium  for  a 
Level-Slice  Process 


The  level  slice  function  is  a nonmonotonic 
nonlinearity  and  hence  we  have  to  use  a nonzero  order  to 
obtain  this  function.  Since  it  has  only  one  sign  change  in 
its  slope  the  first  order  should  be  used.  We  want  to  have 


I 0 


V'in’11  - 


for  I.  < I 
in  a 


for  Xa  1 Tin  < Jb 


for  I.  < I.  . 
b — in 


(4.9) 


2 1 • 

For  simplicity  we  choose  c =-i.  The  form  of  this  function 
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I 


r 


is  shown  in  Fig.  4.8.  With  an  ideal  recording  medium  shown' 

I 

in  Fig.  3.3  and  a=0 , b=l,  the  halftone  screen  density  from 


Fig.  3.6  is 


f (x)  = 


log(la/lr) 


log(Ib/lr) 


0 < x <j 


j < x < 1, 


(4.10) 


as  shown  in  Fig.  4.9.  From  Eq.  (4.8)  the  output  intensity 
for  a general  recording  medium  in  the  first  order  is 


^2trx 

Io(Iin,1)  = I Z } g[logIin'f (x)]e  L dx!2  (4.11) 


and  substituting  for  f(x)  from  Eq.  (4.10)  we  have 


r L/2 

Jokin'11  * Ie{J„  9llogIln-log(Ia/Ir)), 


dx  + 


L , ,2 


gtlogIin"log(Ib/Ir) ]e 


(4.12) 


V^n'11  - ^{gUogiin-iogda/ir)] 
-g[logIin-log(Ib/Ir) ] }2. 


(4.13) 


From  this  equation  we  can  say  that  the  output  intensity  is 
obtained  by  shifting  the  recording  medium  characteristic 
curve  by  amounts  log(Ia/Ir)  and  log(Ib/Ir)  (determined  from 
the  level  slice  characteristic) , subtracting  these  shifted 
curves  from  each  others,  and  then  squaring  the  result. 
This  shows  how  a perfect  level  sl^ce  function  is  obtained 
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when  we  have  a binary  recording  medium  with  a 
characteristic  curve  shown  in  Fig.  3.3.  On  the  other  hand, 
for  a nonbinary  recording  medium  with  a characteristic 
curve  like  Fig.  4.1,  the  toe  and  shoulder  of  this  curve 
makes  the  sharp  corners  of  the  level  slice  function  rounded 
and  the  linear  part  of  this  curve  make  the  rise  and  fall  of 
the  level  slice  function  linear  rather  than  abrupt. 
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CHAPTER  5 


PIECEWISE 

LINEAR 

RECORDING 

MEDIUM 

In  the  previous 

chapter 

we  saw 

how 

the  nonbinary 

characteristic  of 

the  recording 

medium 

affects  the 

nonlinear  process. 

We 

conclude 

that 

for  smooth 

nonlinearities  the  linear  part  of  the  recording  medium 

I 

characteristic  is  the  main  source  of  degradation  for  the 
output  of  the  process. 

In  this  chapter,  we  design  halftone  screen  density 
profiles  that  compensate  for  the  nonbinary  characteristic 
of  the  recording  medium.  This  means  that  we  want  to  find 
f(x)  from  Eq.  (3.5)  once  IQ(Iin,k),  the  form  of  g(.),  and  k 
are  specified.  This  is  a nonlinear  integral  - equation; 

J i 

: problem  which  may  not  have  a general  solution.  Batten  and { 
Everett  [5-1]  in  their  study  of  the  control  of  film! 

i 

i 

| characteristics  have  obtained  closed  form  solutions  to  a 
! similar  integral  equation  for  certain  types  of  nonlinear 
functions.  Unfortunately  their  method  cannot  be  applied  to 
our  formulas  for  tfte  important  cases  of  logarithmic  and 
exponential  processes.  This  occurs  because  the  integral  on 

the  right-hand  side  of  Eq.  (3.5)  is  raised  to  a power  of 

j 

two  while  in  their  equation,  they  only  have  a simple 
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integral.  So,  the  solutions  that  they  have  obtained  for 
the  logarithmic  and  exponential  processes  are  not 
applicable  in  our  formulation. 


To  simplify  somewhat  the  above-mentioned  integral 
equation,  we  assume  that  the  recording  medium  has  a 
piecewise  linear  characteristic  as  shown  in  Fig.  4.4.  This 
assumption,  as  mentioned  above,  is  a good  one  for  monotonic 
nonlinearities.  Noting  that  monotonic  nonlinearities  are 
possible  in  the  zero  order,  we  then  have  k=G  in  Eq.  (3.5) 
leading  to  added  simplification.  With  this  procedure  we 
get  an  integral  equation  which  is  easier  to  work  with  at 
the  expense  of  limiting  ourselves  to  monotonic 
nonl inear ities. 

In  the  followings  we  attempt  to  find  f(x)  from  Eq.(3.5) 
with  the  above  assumptions.  The  results  that  we  obtain 
work  only  for  monotonic  nonlinear  functions.  Compensation 
methods  for  other  forms  of  nonlinearities  will  be  discussed 
in  the  next  chapter. 

5.1  Halftone  Screen  Density  Function  for  Monotonic 
Nonlinearities  Using  a Piecewise  Linear  Recording  Medium 

The  piecewise  linear  recording  medium  characteristic 
curve  is  shown  in  Fig.  4.4.  From  this  figure  we  can  write 


g(* ) = a, 


when  logl.  -f(x)  < logl. 


(5.1) 


g<’>  = l'ogTlT7l V)  { (b~a)  UogIin-f(x)]+alogI2-blogI  }, 

c.  1.  ( J • 

when  logI1£logIin~f (x) <logI2 


g ( • ) = b,  when  log  I-  £ log  I.  -f (x) . (5.3) 


Noting  that  we  are  solving  for  monotonic  nonlinearities, 
f(x)  can  also  be  assumed  to  be  a monotonic  increasing  or 
decreasing  function,  depending  on  the  type  of  nonlinearity. 
Here,  we  assume  f(x)  to  be  a monotonical ly  increasing 
function  which  is  applicable  to  the  increasing  logarithmic 
and  exponential  nonlinearities.  If  we  let  the  inverse 
function  of  f(x)  be  represented  by  h(*)»  then  we  have 


y = 


(5.4a) 


e = logd2/I1) 


(5.4b) 


log(Iin/I2)  = y_E* 


(5.4c) 


Then  the  conditions  of  Eq.  (5.1)  through  Eq.  (5.3)  are 
equivalent  to 


when  x > h(y) 


(5.5) 


(a-b) f (x) +by-a (y-e ) ] , 
when  h(y-e)  < x ^ h(y) 


(5.6) 


and 


g(*)  = b,  when  x < h(y-e) 


(5.7) 


I These  equations  specify  the  characteristics  of  a piecewise 
linear  recording  medium.  Now,  rewriting  the  general 
formula  in  the  zero  order  from  Eq.  (3.6)  in  the  form 

L 
0 

and  using  the  above  requirements  for  g(*)  in  it,  we  get 

h(y-e)  . rb(y) 

bdx+e  I [ (a-b)  f (x)+by-a  (y-e)  ]dx 

(y- 

(5.9) 


T * 

V1™-01  - <eJ  9Uogiin-f<*na*>2  (5.8) 


rL  fhty-ej  , rD iy ; 

g [logl.  -f (x) ] dx  = bdx+±  [ (a- 

Jo  in  '0  e J h(y-e) 


fL 

h (y) 


+ | adx 


i 

| which  after  simplification  becomes 


r 
> 0 


gtlogl.  -f  (x)  ]dx  = ( 

in  E y-e 


h (u)du+aL. 


(5.10) 


The  integral  on  the  right-hand  side  of  the  above  equation. 


I 


can  be  simplified  again  as  described  in  Appendix  A.  The 
result  is 


/•h(y)  r y 

f(x)dx  = yh(y)- (y-e)h(y-e)-  I 
J h (y-e ) Jv- 


and  when  used  in  Eq.  (5.10)  gives 


y 

h (u)  du 
y-e 


(5.11) 


(0  9llogIin"£(x,|d*  ■ II 


h (u)du+aL. 


y-e 


(5.12) 
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This  combined  with  Eq.  (5.8)  reduces  to 


ry 

h (u 

< v-r 


)du+aL]  . 


(5.13) 


Taking  the  square  root  of  both  sides  of  Eq.  (5.13)  gives 


o in 


7TJT  = ^ fY  h(u 


) du+aL. 


(5.14) 


The  positive  sign  has  been  chosen  for  the  square  root  since 
both  y and  h(*)  take  on  nonnegative  values.  The  right-hand 
side  of  Eq.  (5.14)  is  a function  of  I^n,  which  in  turn  is  a 
function  of  y through  Eq.  (5.4).  This  dependence  can  be 
shown  through 


2(y>  = +/l0dln.' 


(5.15) 


and  when  used  in  Eq.  (5.14)  results  in 


= ^15.  [ h(u) du+aL. 

E J V-P 


(5.16) 


To  remove  the  integral,  we  can  differentiate  both  sides  of 
this  equation  with  respect  to  y to  get 


z' (y)  = ^ (h(y)-h(y-e) ] 


(5.17) 


This  is  a difference  equation  in  terms  of  h(*).  One  method 
to  solve  such  an  equation  is  to  use  the  Laplace  transform 
[5-2].  Considering  the  definition  of  y in  Eq.  (5.4),  we 
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see  that  tor  a physically  realizable  problem  this  parameter 
is  nonnegative.  So,  we  can  take  the  Laplace  transform  of 
both  sides  of  Eq.  (5.17)  to  obtain 


Z'  (s)  = ^[H(s)-H(s)e_es] 


(5.18) 


where 


Z'  (s)  = s£  [ z ' (y)  ] 


H (s)  = £ [h  (y)  ] . 


(5.19) 


Equation.  (5.17)  can  be  written  as 


H (s)  = 


e . Z' (s) 

b-a  * . -es 
1-e 


(5. 20) 


This  equation  gives  the  Laplace  transform  of  the  function 
h(*).  Before  transforming  back  to  obtain  h(*),  note  that 


e-es I < 1 


(5.21) 


hence 


-es 


i -es  -2es  -3es 
= 1 + e + e + e +. 


(5.22) 


1-e 

which  permits  Eq.  (5.20)  to  be  written  as 


H (s)  = g^tZ*  (s)+Z'  (s)e"£S+.  . . ] . (5.23) 

Now  taking  the  inverse  Laplace  transform  of  both  sides  of 
Eq.  (5.23)  gives 

h(y)  = gr^t2' (y)+z' (y-e)+z' (y-2e)+. . . ] (5.24) 

57 


which  is  easily  seen  to  satisfy  Eq.  (5.16)  except  for  a 


constant.  When  we  include  this  constant,  the  complete 
solution  is 

h(y)  = b§^[z'  (y)+z'  (y-e)  + - • • 1+gr^  (5.25) 

which  expresses  the  inverse  function  of  the  halftone  screen 
density  in  terms  of  the  derivative  of  the  nonlinear 
function  that  we  want  to  generate.  To  get  a solution  from 
this  formula,  the  nonlinear  function  should  be 
differentiable,  a condition  which  is  true  for  the  usual 
smooth  monotonic  nonlinearities  in  practice.  On  the  other 
hand,  the  solution  for  h(*)\ should  be  monotonic  and  hence 
invertible  to  give  the  density  function.  In  the  following 
sections  we  consider  specific  examples  and  find  out  when  a 
satisfactory  h(*)  exists. 

5. 2 Examples 

In  this  section  we  consider  the  power,  logarithmic, 
and  exponential  transformations  separately.  In  each  case, 
the  desired  input-output  relationship  will  be  used  in  the 
results  of  the  previous  section  to  obtain  the  proper 
halftone  screen  density  function. 


5..  2. 


1 Power  Transformation 

For  this  transformation  we 


want 


In(Iin'0)  = 
o in  i . 


A 


(5.26) 
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where  y and  X are  constants  and  I,  is  defined  in  Fig.  4.4. 
Iin 

Replacing  for  — — in  Eq.  (5.26)  from  Eq.  (5.4)  we  get 
1 


1o(Iin'0>  = Y*10 


Xy 


(5.27) 


and  combining  with  Eq.  (5.15)  gives 

1/2  7^ 

z (y)  = +yx/  -10  2 


(5.28) 


From  this 


z' (y)  = K. • 10 


jXy 


(5.29) 


where 


X 1/2 

— — • V • 


= j-y  ' *£nl0. 


Now  we  use  Eq.  (5.29)  in  Eq.  (5.25)  to  obtain 

eK.  i-Xy  i Xy(y-e)  L 

= r-±[10  2 +10  2 +...]+vF^ 


h(y)  = g-juo 


(5.30) 


(5.31) 


eK.^10 


jXy 


b-a 


— -*-X  £ 

[1+e  +e 


-jX*2e 


+...1.  (5.32) 


Noting  that  the  term  inside  the  bracket  is  a geometric 
series  which  sums  to 


1-e 


_|Xe 


we  can  write 
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r 2Xy  -i 

1 T eKle 

h(y)  _ r—  +aL  • 

L "2  Ae  J 


(5.33) 


This  can  easily  be  inverted  to  obtain  f(x)  as 


CM  - 2 2 )). 


(5.34) 


5.2.2  Logarithmic  Transformation 


For  this  transformation  we  want 


ioUin-O)  - K-logd^/Ij) 


(5.35) 


where  K is  a constant  and  I is  defined  in  Fig.  4.4. 

Iin 

Replacing  for  — in  Eq.  (5.35)  from  Eq.  (4.4)  we  get 
1 1 


I0(Iin.°>  - Ky 


(5.36) 


and  combining  with  Eq.  (5.15)  gives 


* (y)  = +(Ky)1/2. 


(5.37) 


From  this  we  obtain 


.•(»>-  y/2y-1/2. 


When  used  in  Eq.  (5.25)  this  gives  us 


(5.38) 


h(y)  " iy~1/2+<y-e>~1/2+- 


(5-39) 


Noting  that  z'(y)  is  a monotonically  decreasing  function  we 
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see  that  h(y)  is  not  necessarily  monotonic.  Hence,  in  this 
case,  we  cannot  invert  h(y)  to  obtain  f(x),  although  it 
will  be  obtained  in  the  next  chapter  by  another  method. 

5.2.3  Exponential  Transformation 


For  this  transformation  we  want 

I . 
in 


loam.0'  - «<B>  *1 


(5.40) 


where  a and  8 are  constants  and  I-,  is  defined  in  Fig.  4.4, 
I . 

Replacing  for  -j — in  Eq.  ( 5 . 40 ) from  Eq.  (5.4)  we  have 
1 


V^n'01  = a(B) 


10J 


(5.41) 


and  combining  with  Eq.(5.15)  gives 


z (y)  « +a1//2  (6) 1//2’ 1Q} 


(5.42) 


From  Eq.  (5.42)  we  have 


z’ (y)  = K9*10y*B10 


(5.43) 


where 


K2  = J a^2*  £nl0*  £nB. 


(5.44) 


Now  we  use  above  in  Eq.  (5.25)  to  obtain 


h(y) 


j~10y[eloY+l<fe' 


,ioy-e+. 


1 + 


aL 

b-a. 


(5.45) 
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It  can  be  seen  that  z'(y)  is  a monotonical ly  increasing 
function  which  would  make  h(y)  a monotonically  increasing 
function  also.  Although  we  cannot  invert  h(y) 
analytically,  it  is  possible  to  do  so  by  numerical  means. 

5. 3 Comments 

Equation  (5.25)  gives  the  inverse  function  of  the 
halftone  screen  density  function  for  a piecewise  linear 
recording  medium.  To  obtain  the  density  function  we  should 
be  able  to  invert  the  function  obtained  from  this  equation. 
The  logarithmic  process  example  shows  that  this  inversion 
is  not  possible  in  general.  This,  of  course  does  not  mean 
that  the  compensation  problem  for  monotonic  nonlinearities 
should  be  restricted  to  a certain  class  of  such  functions. 
The  fact  is  that  too  many  restrictions  are  placed  on  our 
problem  to  make  it  manageable  with  the  mathematical 
formulation  presented  here.  In  the  next  chapter  we 
consider  the  compensation  problem  in  general  and  try  to 
find  new  methods  of  compensating  for  the  nonbinary 
characteristic  of  the  recording  medium  for  both  monotonic 
and  nonmonotonic  nonlinearities. 
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CHAPTER  6 


DISCRETE  DENSITY  HALFTONE  SCREEN 

Several  methods  are  available  in  practice  for 
generating  halftone  screens  with  a desired  density  profile. 
Some  methods  are  purely  optical  and  involve  the 
photographic  recording  of  geometrical  shadows  or 
diffraction  patterns  from  ruled  gratings.  Although  this 
technique  produces  continuous  halftone  screens,  it  does  not 
have  the  flexibiJ ity  to  precisely  produce  arbitrary  screens 
needed  for  nonlinear  processing. 

Another  type  of  halftone  screen  density  profile  that 
can  be  generated  in  practice  is  a step  function 
approximation  to  the  desired  continuous  density.  These 
halftone  screens  are  generated  by  digital  image  recorders, 
plotting  microdensitometers  or  step  and  repeat  cameras. 
Hence  the  theoretical  accuracy  available  in  designing  the 
halftone  screen  density  profile  is  limited  by  the  practical 
limitations  in  making  the  screen.  This  motivates  the 
following  analysis  which  considers  the  halftone  screen 
density  profile  as  a step  function  approximation  with 
discrete  values.  As  will  be  shown  in  this  chapter,  this 
assumption  helps  to  simplify  the  formulas  and  allows  us  to 
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obtain  some  conclusions  which  cannot  be  drawn  from  the 
continuous  density  analysis.  To  do  so,  we  formulate  with  a 
discrete  density  halftone  screen  and  derive  optimization 
formulas  for  the  zero  and  nonzero  orders. 


6, 1 Zero  Order 

The  general  formula  relating  the  output  intensity  in 
the  zero  order  to  the  input  intensity  is  given  by 

rL 


Io(Iin'0)  = {e  (ogn°giin-f(x)]dxr 


(6.1) 


as  expressed  in  Eq.  (3.6).  Let  f(x)  be  approximated  as 
shown  in  Fig.  6.1,  then  Eq.  (6.1)  can  be  written  as 

Io(Iin'0)  = 1g(logIin-a1)dx+j  2 g (logIin-a2)dx+ 

j 0 J L, 

(6.2) 


■A 


g(logIin~aN)dx] } * 


JN-1 


Assuming  that  we  have  equal  dimensions  on  the  x-axis,  i.e, 


Ll_0  = L2-Ll  * ‘ ' LN  LN-1  N 


(6.3) 


then 


Io(Iin'0)  = {E[ (Ll-0)g(logIin'al)+(L2_Ll)g(logIin“12)+ 

(6.4) 

. . .+(LN-LN_1)g(logIin~aN) ] } 


or 
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obtain  some  conclusions  which  cannot  be  drawn  from  the 
continuous  density  analysis.  To  do  so,  we  formulate  with  a 
discrete  density  halftone  screen  and  derive  optimization 
formulas  for  the  zero  and  nonzero  orders. 


6-.  1 Zero  Order 

The  general  formula  relating  the  output  intensity  in 
the  zero  order  to  the  input  intensity  is  given  by 


I (I -n,0)  = (i 
o in  u 


g [logIin~f (x) ] dx} ' 


(6.1) 


as  expressed  in  Eq.  (3.6).  Let  f(x)  be  approximated  as 
shown  in  Fig.  6.1,  then  Eq.  (6.1)  can  be  written  as 


Id.,0)  = M 

o m i* 


ri 


g dogIin"a1)  dx+ 


g (logl . -addx+ 
111  ^ 


(6.2) 


g(logIin"aN)dxl } 


JN-1 


Assuming  that  we  have  equal  dimensions  on  the  x-axis,  i.e. 


Ll-0  - L2-Ll  = * * ‘ = LN  LN-1  N 


(6.3) 


then 


1o(1in'0)  = ^l(Ll'0)g(logIin"al)  + (L2_Ll)g(logIin_a2)  + 

(6.4) 

••'+(LN-LN-l)g(logIin"aN)1}2 


or 
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N 

Io(Ii„'0)  * {E  Zg(logIln-aN))2.  (6.5) 

The  above  formula  gives  the  output  intensity  in  the  zero 
order  as  a function  of  the  discrete  grey  levels  on  the 
halftone  screen  and  the  characteristic  curve  of  the 
recording  medium. 

6.1.1  Synthesis  Problem  (Zero  Order) 

The  synthesis  problem  here  is  to  find  a i=l,N  in 

Eq.  (6.5)  when  g(»)  and  the  functional  relationship  of 

Io^Iin'tJ^  *n  terms  of  *s  9iven*  Since  Eq.  (6.5)  holds 

over  the  input  intensity  range,  the  functional  relationship 

of  in  terms  of  I^n  holds  for  discrete  values  of 

I.  , say  I.  =1.  , where  k=l,K.  Thus  in  discrete  form 
in  in  it 

Eq.  (6.5)  can  be  written  as 

N 

W01  ■ 
j N 

Io(Ik'0)  = Sgdogik-ai)}2  (6.6) 

: T 

One  procedure  to  find  the  optimum  a^ ’ s is  to  minimize  the 
mean  square  error  expression 

i ...  - 66  j 

I 

. 

/ 


gl . -a.T 
^ in  N 


The  above  formula  gives  the  output  intensity  in  the  zero 
order  as  a function  of  the  discrete  grey  levels  on  the 
halftone  screen  and  the  characteristic  curve  of  the 
recording  medium. 


6.1.1  Synthesis  Problem  (Zero  Order) 


The  synthesis  problem  here  is  to  find  a^,  i=l,N  in 
Eq.  (6.5)  when  g(*)  and  the  functional  relationship  of 
Io(I^n,0)  in  terms  of  Iin  is  given.  Since  Eq.  (6.5)  holds 
over  the  input  intensity  range,  the  functional  relationship 
of  terms  of  I^n  holds  for  discrete  values  of 
If  , say  where  k=l,K.  Thus  in  discrete  form 
Eq.  (6.5)  can  be  written  as 


IQ(I1,0)  - (jy  Sg(logI1-a1)): 

j T 

louk,0>  = Eg(logIk-ai)}: 

: ^ 

■o1^'01  * Lgdogin-ai))1 


(6.6) 


One  procedure  to  find  the  optimum  a^'s  is  to  minimize  the 
mean  square  error  expression 


min  E 


K N 

= 2 {I  (l.,0)-[i  2g(logI  -a.)]2}2.  (6.7) 
k=l  ° i=l  K x 

ai'a2 ' • • • 'aN 

| i 

This  should  produce  values  of  ads  that  bring  the  output 
intensity  in  the  zero  order  as  close  as  possible  to  the 

desired  output  intensity  in  the  mean-square  sense. 

I 

Now,  since  a^,...,aN  are  the  different  density  values 

i 

on  the  halftone  screen,  their  values  should  be  positive  and 
remain  between  certain  practical  limits,  as  expressed  by 

i 

) 

c 1 ai  1 d (6.8) 

i 

where  c and  d are  given  constants.  This  makes  the 

minimization  problem  a constrained  minimization.  There  is 

j 

a transformation  that  simplifies  this  problem  and  imbeds 

I 

the  constraint  in  the  expression  to  be  minimized  [6-1]  . If 
we  let 

i 

l 

a^  = (d-c) sin2y^+c  (6.9) 

this  makes  the  values  of  aiSlimited  to  the  range  [c,d]. 
Note  that 

3F  _ 3E  8ai  i - i m 
^-^7^y7  1 " 1'N-  (6.10) 

When  we  minimize  E with  respect  to  y , we  set 

- 0 i « 1,N  (6.11) 

L. 


67 


and  combining  Eq.  (6.10)  and  (6.11)  we  have 


or 


(6.12) 


3a. 

— - =0  i = 1,N.  (6.13) 

3yi 

To  perform  the  minimization  we  want  Eq.  (6.12)  to  be 
satisfied  at  the  same  time  that  Eq.  (6.13)  is  not  satisfied 
in  the  range  [c,d].  To  check  for  this  condition  note  that 

33i 

5y-  = (d-c)sin2yi  (6.14) 

which  when  set  to  zero  gives 

yi  = 0 =>  = c 

or  i = 1,N  (6.15) 


It  can  be  seen  that  the  only  values  of  a!s  in  the  range 

3a^  1 

[c,d]  that  make -gy—  =0  are  the  boundaries  of  this  interval, 
and  this  can  be  prevented  from  causing  any  problem  by 
choosing  c and  d,  such  that  the  interval  [c,d]  contains  the 
limit  values  for  the  a^'s. 

6.1.2  Example  of  A Function  Possible  with  the  Zero  Order: 
Logarithmic  Screen 

In  a logarithmic  process  we  want  the  relationship 


Io(Iin'0)  = ml°g(Iin/Ir)  (6*16) 
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min  E = £ tIo(Ik'0)_[^  2g(logIk-ai)  ]2}2.  (6.7) 
k—1  i“ 1 

ai'a2 ' • • • 'aN 

This  should  produce  values  of  ai’s  that  bring  the  output 
intensity  in  the  zero  order  as  close  as  possible  to  the 

i 

desired  output  intensity  in  the  mean-square  sense. 
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i 

c < a.  < d (6.8) 

— X 

where  c and  d are  given  constants.  This  makes  the 

minimization  problem  a constrained  minimization.  There  is 

l 

a transformation  that  simplifies  this  problem  and  imbeds 

the  constraint  in  the  expression  to  be  minimized  [6-1],  If 

we  let 

| j 

a^  = (d-c) sin2y^+c  (6.9) 

this  makes  the  values  of  aislimited  to  the  range  [c,d]. 
Note  that 

3E  _ 3E  9ai  , _ , N 

3y7  ~ -5¥7  ?y7  1 " 1'N‘  (6.10) 

When  we  minimize  E with  respect  to  y , we  set 

#-  = 0 i - 1,H  (6.11) 

9yi 


and  combining  Eq.  (6.10)  and  (6.11)  we  have 


or 


_9E_ 

9ai 


3Yi 


= 0 


i = 1 , N 


(6.12) 


= 0 


i = 1 , N. 


(6.13) 


To  perform  the  minimization  we  want  Eq.  (6.12)  to  be 
satisfied  at  the  same  time  that  Eq.  (6.13)  is  not  satisfied 
in  the  range  [c,d].  To  check  for  this  condition  note  that 


3a . 


— = (d-c) sin2y . 


3y 

-L 

which  when  sec  to  zero  gives 

Yi  = 0 =>  ai  = c 

TT  . 

Yi  = _ =>  a.  = d 


(6.14) 


or 


i = 1 , N 


(6.15) 


It  can  be  seen  that  the  only  values  of  a|s  in  the  range 

3a^ 

[c,d]  that  make-gy7=0  are  the  boundaries  of  this  interval, 
and  this  can  be  prevented  from  causing  any  problem  by 
choosing  c and  d,  such  that  the  interval  [c,d]  contains  the 
limit  values  for  the  a^'s. 

6.1.2  Example  of  A Function  Possible  with  the  Zero  Order: 
Logarithmic  Screen 

In  a logarithmic  process  we  want  the  relationship 


I<d  ( Jin ' 0)  = 


(6.16) 
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I 


between  the  output  and  input  intensities  where  m is  a 


r 


i 

/ 

i 


I 


constant  and  Ir  is  shown  in  Fig.  3.3.  The  corresponding 
screen  density  profile  for  a binary  recording  medium  shown 
in  Fig.  3.3  is  obtained  from  Eq.  (3.19)  as 

f (x)  = i[a+^-x]2,  0 < x < L.  (6.17) 

Hence  the  initial  values  for  the  a.'s  are 

i 

a.  = — [a+^jr^  x.]2,  0 < x.  < L,  i = 1,N.  (6,18) 

l m L i — l — 

From  Eq.  (6.9)  the  corresponding  y^ ' s are 


Y±  = sin 


-1/1  1 r b-a  , 2 

sIa+— -° 


(6.19) 


Hence  to  find  the  optimum  screen  giving  a logarithmic 
relationship  as  in  Eq.  (6.16)  with  a recording  medium 
characteristic  curve  g(logE),  we  must  perform  the 
minimization 


min  E 


yl,y2 


»y 


N 


^2  {mlog(Ik/Ir)-[i]Cg(logIk-(d-c) 
k=l  i=l 

2 n2,2 

sin  yi-c) J ) 


(6.20) 


with  the  initial  values  of  Eq.  (6.19). 

A computer  routine  has  been  generated  which  performs 
the  above  minimization.  A subroutine  which  is  taken  from 
the  IMS L library  [6-2]  and  uses  a quasi-Newton  algorithm 
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for  finding  the  minimum  of  a function  of  N variables  [6-3] 
is  the  main  part  of  this  routine.  The  results  for  a 

! 

piecewise  linear  recording  medium  (Fig.  4.4  with  a=0,  b=l) 

: . I 

with  different  slopes  in  its  linear  part  (gamma)  are  shown| 
in  Figs.  6.2  through  6.9.  It  is  assumed  that  there  are  30 

t 

discrete  points  in  the  halftone  screen  density  profile  and 

I 

I 

the  density  values  are  between  0 and  2 in  these  figures. 
It  is  also  assumed  that  the  values  of  log  I (shown  in 

i r 

Fig.  3.3)  lie  in  the  middle  of  the  values  of  log  ,and  log 

i . 

1^  (shown  in  Fig.  4.4).  Note  that  the  plots  of  the 

input-output  curves  are  semi-logarithmic  and  hence  the 

ideal  result  is  a straight  line.  The  optimized  output 

curve  is  seen  to  approximate  the  ideal  result  with  a small 

error . 

To  examine  how  the  discretization  affects  the 
degradation  of  the  output,  we  have  computed  the  degradation 
from  Eq.  (6.5)  for  a recording  medium  as  above  with  a gamma 
equal  to  one  and  for  different  numbers  of  points  in  the 
halftone  screen  density  profile.  The  result  is  shown  in 
Fig.  6.10.  Figure  4.6b  is  a similar  result  with  a 
continuous  density  profile,  and  comparing  it  with  Fig.  6.10 
shows  that  for  higher  values  of  N the  two  graphs  agree 
quite  well. 

6. 2 Nonzero  Order 

The  general  formula  relating  the  output  intensity  in 
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between  the  output  and  input  intensities  where  m is  a 
constant  and  Ir  is  shown  in  Fig.  3.3.  The  corresponding 
screen  density  profile  for  a binary  recording  medium  shown 
in  Fig.  3.3  is  obtained  from  Eq.  (3.19)  as 

f(x)  = i[a+^x)2,  0 < x < L.  (6.17) 

Hence  the  initial  values  for  the  a.'s  are 

i 

a.  = i[a+^  x.]2,  0 < x.  < L,  i = 1,N.  (6.18) 

-L  m l 1 — l — 

From  Eq.  (6.9)  the  corresponding  y^ ' s are 


yi  " sin  Va^E  Sla+TTxi]2-c  • 


(6.19) 


Hence  to  find  the  optimum  screen  giving  a logarithmic 
relationship  as  in  Eq.  (6.16)  with  a recording  medium 
characteristic  curve  g(logE),  we  must  perform  the 
minimization 


min  E 


y,  >y 


2' * * * ,yN 


N N 

{mlog(Ik/Ir)-[iX]g(loglk-(d-c) 

k=l  i=i 

2 .,2,2 

sin  y^c)  J } 


(6.20) 


with  the  initial  values  of  Eq.  (6.19). 


A computer  routine  has  been  generated  which  performs 
the  above  minimization.  A subroutine  which  is  taken  from 
the  IMSL  library  [6-2]  and  uses  a quasi-Newton  algorithm 
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! 

discrete  points  in  the  halftone  screen  density  profile  and 
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the  density  values  are  between  0 and  2 in  these  figures. 

It  is  also  assumed  that  the  values  of  log  I (shown  in 
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i 

input-output  curves  are  semi-logarithmic  and  hence  the 
ideal  result  is  a straight  line.  The  optimized  output 
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| 

To  examine  how  the  discretization  affects  the 
degradation  of  the  output,  we  have  computed  the  degradation 
from  Eq.  (6.5)  for  a recording  medium  as  above  with  a gamma 
equal  to  one  and  for  different  numbers  of  points  in  the 
halftone  screen  density  profile.  The  result  is  shown  in 
Fig.  6.10.  Figure  4.6b  is  a similar  result  with  a 
continuous  density  profile,  and  comparing  it  with  Fig.  6.10 
shows  that  for  higher  values  of  N the  two  graphs  agree 
quite  well. 

6. 2 Nonzero  Order 

The  general  formula  relating  the  output  intensity  in 
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Figure  6.2  Logarithmic  transfer  function  for  a piece- 
wise  linear  recording  medium  with  gamma=1.0 

(a)  ideal;  (b)  degraded;  (c)  optimized 


f(x) 


Figure  6.9  Halftone  cell  shape  corresponding  to  Figure  6.8 
(a)  ideal;  (b)  optimized 
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Figure  6.10  Logarithmic  transfer  function  for  a piecewise 
linear  recording  medium  with  gamma=1.0 

(a)  ideal 

(b)  16  points  in  the  discrete  density  profile 

(c)  200  points  in  the  discrete  density 
profile 
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the  nonzero  orders  to  the  input  intensity  is 


L j2?rkx 

[ (I.  ,k)  =|{fL  gtlogl.  -f(x)]e  Ldx}|2  (6.26) 

u in  ij  q in 


from  Eq.  (3.5a).  Using  a quantized  approximation  to  f(x), 
as  in  the  zero  order  case  we  can  write  (from  Fig.  6.1) 


i A Li 


. 2irkx 


Io(Iin'k)  = Il[  0 5(1°glin-a1)e  dx+ 


(6.27) 


.2irkx  T .2nkx 

fL  ^ Ij  - 2 

g(logI.  -a9)e  dx+...+  g(logI.  -a  )e  dx] | 

11*1  Cm  _ in  n 

L1  ln-i 

For  simplicity  we  let 


g(logIin-ai)  = g.^  i = 1,N. 


(6.28) 


Then 


VW10 


. 1,2  I gi 

4 tt  k 


2irkL^  2nk*  0 


2itkL 


11  L 3-T—  + . ] — 

-e  +. . . +gM  (e 

j’-U  2 

e L ) |2 


~T71  lgl* 

4,t  k 


■ 2irk«  0 .27Tk(Ll  0) 
e L (e3  L -1)+...+ 


2nkLN_1  2Trk(L-LN_1) 


(6.29) 


1 f 7 

gNe  (e 


Now  from  Eq.  (6.3) 


Ll-0  = L2~L1  = L3~L2  = • 


JN-1  N' 


(6.30) 


so  when  Eq.  (6.30)  is  used  in  Eq.  (6.29)  we  have 
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2iTkL. 


2-rrkL, 


o in 


4tt  k 


g1+g2e 


+...+gNe 


Note  also  that 


i N -I  | 2 

|e  -1|  . 


(6.31) 


T,  _ L t — 9— 

L1  ~ N'  2 ZN'  ' 


* Li  - i 


(6.32) 


consequently  Eq.  (5.31)  can  be  written  as 


I (I.  , k)  = — i 
° in  47r2k2 


N-1|2-I  £g(iogiin-a.)e 


. 2irk  ( i-1 ) 
] L 


(6.33) 


Equation  (6.33)  gives  the  output  intensity  in  any  nonzero 
order  k in  terms  of  the  input  intensity,  the  characteristic 
function  of  the  recording  medium  and  the  discrete  grey 
levels  on  the  halftone  screen. 


6.2.1  Synthesis  Problem  (Nonzero  Order) 


The  synthesis  problem  here,  as  in  the  case  of  the  zero 
order,  is  to  determine  suitable  ai  values  from  Eq.  (6.33). 
The  procedure  that  we  take  is  the  same  as  the  one  for  the 
zero  order.  Namely,  we  minimize  the  expression 


min  E 

al'a2' * ' ’ ,aN 


M j2nk 

- £ ,VIin-k)-r47le  " -1'2- 

m=l  4 it  k 

N .2Trk(i-l) 

I 2 g(logIm-ai)e  L |2]2. 


(6.34) 
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To  limit  the  resulting  density  values  as  in  Eq.  (6.8)  we 
transform  the  problem  to  the  minimization 


min  E 

yl'y2' * * ’ ,yN 


= £ (I  (I,  ,k) T~l 

m=l  ° ln  4irV 


. 2uk 

3 L . | 2 
e -1  . 


(6.35) 

N j2irk  (i-1) 

| X]g [logl^- (d-c) sin2yi-c] e L |2}2. 


The  yi,s  are  related  to  the  a i’ s through  Eq.  (6.9).  To 
initialize  the  values  of  the  y s for  this  minimization 
procedure,  we  use  the  density  values  obtained  for  the 
screen  when  the  recording  medium  is  binary. 

6.2.2  Example  of  a Function  Possible  in  the  Nonzero  Order: 
Level  Slice  Screen 

To  obtain  a level  slice  transformation  the  first  order 
is  a suitable  choice.  For  this  function  we  want 


V'in'11  * 


for  I.  <1, 
in  a 

for  I < I.  < IK 
3 — in  b 

for  I.  > I. 
in  — b 


(6.36) 


and  this  function  is  shown  in  Fig.  4.6.  The  density 
profile  giving  such  a relationship  for  the  binary  recording 
medium  shown  in  Fig.  3.3  is 


f(x)  = 


log ( Ia/Ir)  for  0 < x < j 


log(Ib/Ir)  for  j<  x< 


(6.37) 


from  Eq.  (4.10).  The  corresponding  density  levels  in  the 
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I 


L 


v 


discrete  screen  are 

log  (I  /I  ) for 
d r 


ai  = 


• ^ N 
1 < 2 


N 


log(Ib/Ir)  for  i > ^ • 
Using  Eq.  (6.38)  in  Eq.  (6.9),  we  have 


(6.38) 


■ 


. -i 
sm  " 


-c 


d-c 


_ . N 

for  1 < 2 


\ sm 


- 2 


(6.39) 


, . N 

1 ^r  1 > o - 


With  the  initial  values  for  y^ ' s as  above  we  then  want  to 
minimize 


M 


2ir 

3 M 


min  E =2  (I  (I  ,D — iy|eJN-l[2 


, o m . « 

m=l  4 7r 


N 


i=l 


(6. 40) 


g [logIm~  (d-c)  sin2y^-c]  | 2 } 2 


A computer  routine  similar  to  the  one  for  the 
logarithmic  process  has  been  generated  to  perform  the  above 
minimization.  With  the  same  assumptions  about  the 
recording  medium  and  the  halftone  screen  as  with  the 
logarithmic  process  (section  6.1.2),  the  result  is  shown  in 
Figs.  6.11  through  6.18.  It  is  seen  that  the  optimization 
procedure  is  successful  in  restoring  the  shape  of  the  level 
slice  function,  but  cannot  do  much  for  the  edges  with 
finite  slope  in  the  degraded  response.  As  explained  in 
section  (4.2.1),  the  finite  slope  edges  are  a direct 
consequence  of  the  characteristic  curve  of  the  recording 
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CHAPTER  7 

EXPERIMENTS  WITH  HALFTONE  SCREENS 

7. 1 Halftone  Screen  Fabrication 

To  generate  a halftone  screen  we  must  plot  a given 
density  profile  (such  as  those  obtained  in  Chapter  6)  on 
photographic  film.  To  do  so,  various  methods  have  been 
tried  in  previous  work.  Some  more  recent  methods  are: 
optical  filtering  of  binary  gratings  by  Strand  (7-1],  use 
of  a plotting  microdensitometer  by  Dashiell  and  Sawchuk 
[2-1],  and  moving  a slit  grating  on  top  of  the  photographic 
film  with  repeated  exposures  as  done  by  Liu,  Goodman  and 
Chan  [7-2] . 

We  have  used  a computer  controlled  Dicomed  Digital 
Image  Recorder  to  generate  our  screens.  This  machine  can 
plot  up  to  4096x4096  discrete  points  in  a 4"x4"  area. 
There  are  256  different  intensity  levels  for  each  exposure 
setting.  To  plot  any  pattern  on  the  film  with  this 
machine,  an  input  matrix  of  numbers  between  zero  and  255 
representing  the  digitized  pattern  is  read  into  the 
computer.  Before  doing  so,  we  have  to  find  out  what  is  the 
relationship  between  those  numbers  and  the  final  resulting 
density  on  the  film.  To  find  out  this  relationship  a grey 


85 


level  test  pattern  was  plotted  by  the  Dicomed  on  Kodak 
SO-115  film  (the  same  film  used  for  the  screens) . The 
exposed  film  was  then  developed  for  the  recommended  time  in 
D-19  developer.  The  densities  on  the  test  pattern  were 
then  measured  on  a Macbeth  transmission  densitometer.  This 
procedure  was  repeated  for  different  exposure  settings  on 
the  Dicomed.  That  results  in  a series  of  curves  which 
represents  the  final  density  on  SO-115  film  versus  the 
input  numbers  to  the  computer  for  each  exposure  setting  on 
the  Dicomed.  From  these  curves  one  that  covered  the 
density  range  of  our  halftone  screen  was  selected.  This 
curve  was  then  used  as  the  basis  for  fabrication  of  that 
screen. 

As  a check  on  the  calibration,  a test  pattern  was 
generarated  on  a separate  film  with  the  same  exposure 
setting  of  the  Dicomed  each  time  a new  screen  was  made.  If 
any  deviations  were  observed,  the  exposure  setting  was 
re-adjustcd  and  the  experiment  was  repeated. 

7.2-  Experiments  with  the  Logarithmic  Screen 

To  make  a logarithmic  screen  experimental  data  for  the 
characteristic  curve  of  the  standard  LCLV  (45°  twisted 
nematic  hybrid  field  effect)  was  gathered.  This  data  was 
then  used  in  the  optimization  routines  for  the  logarithmic 
process  (sec.  6.1.2)  to  obtain  the  density  profile  of  the 

optimized  screen.  The  corresponding  halftone  screen  was 
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later  generated  using  the  method  described  in  section  7.1. 
The  screen  was  then  used  in  a real-time  experiment  to  test 
its  effectiveness. 

7.2.1  Logarithmic  Screen  Fabrication 

The  experimental  set  up  to  determine  the  LCLV 
characteristic  curve  is  shown  in  Fig.  7.1.  With  no  input, 
the  control  surface  of  the  light  valve  is  illuminated 
directly  by  an  incoherent  arc  lamp  source.  Readout  was 
performed  coherently  at  632.8  nm  with  a simple  optical 
filering  system.  In  this  case,  no  spatial  filter  was  used 
in  the  filter  plane.  The  resulting  output  intensity 
transmittance  versus  input  intentisy  is  shown  in  Fig.  7.2. 
By  taking  the  square  root  of  the  output  intensity 
transmittance  valves,  and  the  logarithm  of  the  input 
intensity  values  a plot  of  the  amplitude  transmittance  vs. 
log  exposure  curve  of  the  device  was  obtained.  This  is 
shown  in  Fig.  7.3  after  normalization. 

A discrete  version  of  this  graph  was  then  fed  to  the 
computer  routine  for  generation  of  the  optimized  density 
profile.  The  result  is  shown  in  Fig.  7.4.  There  are  16 
discrete  points  in  this  graph.  Using  these  values  in  the 
procedure  of  section  7.1,  the  corresponding  halftone  screen 
was  made.  Recalling  that  we  can  plot  4096  points  over  86 
mm  with  this  method,  the  halftone  screen  has 
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Figure  7.1  Experimental  set-up  for  real-time  parallel  logarithmic 
filtering 


x b4 — = 3 lines/mm 
16  86  iron 

This  halftone  screen  was  then  placed  in  the  input 
plane  of  the  experimental  setup  shown  in  Fig.  7.1.  Note 
that  the  logarithmic  screen  was  designed  on  the  assumption 
of  using  the  zero-order  of  the  diffraction  pattern  of  the 
halftoned  picture.  A simple  pinhole  spatial  filter  at  the 
position  of  the  zero  order  in  the  filter  plane  was  used  to 
pick  out  this  order.  The  resulting  transfer  function  is 
shown  in  Fig.  7.5.  It  approximates  a logarithmic  transfer 
function  with  less  than  5%  error  over  one  decade  and  less 
than  10%  error  over  another  decade. 


7.2.2  Logarithmic  Screen  Performance 


To  test  the  effectiveness  of  the  logarithmic  screens 
in  another  experiment,  two  crossed  multipl icatively 
combined  Ronchi  rulings  were  used  as  an  input  picture  for 
the  experimental  setup  of  Fig.  7.1.  The  period  of  these 
rulings  was  approximately  3 mm,  much  higher  than  the 
halftone  screen  period  of  0.33  mm.  The  spectrum  in  the 
filter  plane  is  shown  in  Fig.  7.6a.  Next,  the  logarithmic 
halftone  screen  was  placed  in  contact  with  the  rulings. 
The  filter  plane  spectrum  is  shown  in  Fig.  7.6b.  The 
difference  in  Fourier  spectra  between  multipl icatively 
combined  gratings  and  additively  combined  gratings  obtained 
by  real-time  logarithmic  filering  is  as  follows.  The 
additive  spectrum  components  lie  only  on  the  x and  y axis 


Figure  7.5  Transfer  characteristic  of  the  LCLV 
with  the  optimized  screen 
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Oth 


2nd 


Figure  7.6  Fourier  spectra  of  crossed  gratings 


(a)  multiplicative  spectrum 

(b)  additive  spectrum  (result  in  the 
Oth  order) 
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around  the  zeroth  diffraction  order  in  the  frequency  domain 
while  the  multiplicative  spectrum  contains  cross-term 
off-axis  components.  Figure  7.6  also  shows  the  higher 
diffraction  orders  that  arise  due  to  halftone  screen.  For 
simple  logarithmic  processing,  these  higher  orders  would  be 
eliminated  by  spatial  filtering. 

The  two-dimensional  point  logarithm  obtained  here  in 
real  time  is  particularly  useful  in  a signal  processing 
technique  called  homomorphic  filtering  [7-3, 7-4]. 
Homomorphic  filtering  is  useful  for  processing  multiplied 
or  convolved  signals,  with  specific  applications  in 
radiography  [7-4, 7-5],  speech  processing,  and  image 
processing  [7-3] . 


CHAPTER  8 


DIRECT  METHOD  FOR  OBTAINING  NONLINEAR  FUNCTIONS 

IN  REAL  TIME 

». 

Although  the  halftone  screen  method  is  a flexible  way 
to  obtain  different  nonlinear  functions,  it  is  often 
desirable  to  bypass  the  halftoning  step  and  obtain  the 
nonlinearities  directly.  By  eliminating  the  halftone 
screen  we  rely  on  the  nonlinear  characteristic  of  the 
recording  medium  to  generate  an  overall  nonlinearity.  Tai, 
Cheng  and  Yu  [8-1]  have  used  the  nonlinear  portion  of  the 
photographic  film  characteristic  curve  to  generate  a 
logarithmic  function.  Because  their  method  uses 
photographic  film,  it  cannot  be  implemented  in  real  time. 

In  this  chapter,  we  discuss  how  the  nonlinear 
characteristic  of  a special  birefringent  liquid  crystal 
device  can  be  used  to  produce  nonlinear  transformations  in 
real  time.  In  particular,  we  consider  implementation  of 
the  analog-to-dig ital  conversion  nonlinearity  by  this 
method . 
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6 • Biref r inqent  Liquid  Crystal  Device 

The  liquid  crystal  light  valve  described  in  chapter 
two  can  be  used  as  a pure  birefringent  device.  In  this 
case  no  twist  is  present  in  the  alignment  of  the  liquid 
crystal  with  its  substrates.  As  the  input  light  changes 
the  birefringence  of  the  liquid  crystal  layer  through  the 
photoconductor,  the  output  beam  goes  through  a series  of 
phase  changes.  This,  in  turn,  is  transformed  into  output 
intensity  variations  by  the  crossed  polarizers.  The  output 
intensity  variaion  with  respect  to  the  applied  voltage  to 
the  liquid  crystal  (no  input  light)  is  shown  in  Fig.  8.1, 
and  its  variation  with  respect  to  the  input  intensity  is 
shown  in  Fig.  8.2.  It  is  seen  that  the  optically 
controlled  device  characteristic  curve  has  fewer  intensity 
peaks  than  the  electrically  controlled  one.  This  is  due  to 
the  saturation  of  the  photoconductor  material  at  high  input 
light  intensities. 

The  direct  nonmonotonic  nature  of  the  device  is  useful 
for  some  limited  approximations  to  obtain  nonmonotonic 
functions.  Different  parts  of  the  curve  shown  in  Fig.  8.2 
can  approximate  different  polynomials.  Each  period  of  this 
curve  also  approximates  a sine-squared  function.  A 
particularly  useful  function  is  the  optical  parallel 
analog- tc-dig ital  conversion  function  which  is  discussed  in 
more  detail  in  the  following. 
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8.2  Real-Time  Parallel  Analog-to-Diqital  Conversion 


The  process  of  analog-to-dig ital  (A/D)  conversion  is 
the  representation  of  analog  (continuous)  information  in 

i 

digital  (binary)  form.  For  one-dimensional  signals  any 
coding  procedure  which  assigns  to  each  signal  value  a group 
of  digits  (bits)  performs  such  an  operation.  For 

! 

j two-dimensional  signals  (images)  we  must  find  a method 
which  assigns  to  every  continuous-tone  image  a number  of 

I 

binary  images  (bit  planes) . Each  binary  image  records  the 
value  of  one  particular  bit  of  the  quantized  intensity  at 
each  image  point.  Fig.  8.3  shows  how  the  A/D  conversion 
j bit  planes  may  be  obtained  from  the  device  characteristic. 
For  simplicity,  it  is  assumed  that  the  device  has  a uniform 
sine  squared  characteristic  over  two  periods  as  shown  in 
Fig.  8.3a.  By  thresholding  the  output  at  1/2  we  get  a 

binary  function  wich  can  serve  as  one  bit  of  a three  bit 

I 

binary  code.  The  other  two  bits  are  then  obtained  by i 

I i 

attenuating  the  input  intensity  to  effectively  get  the! 
first  cycle  and  first  half-cycle  of  the  characteristic 
curve.  These  curves  are  then  thresholded  at  the  1/2  level 
as  shown  in  Figs.  8.3b  and  c.  Note  that  the  binary 
functions  in  Fig.  8.3  are  the  three  bits  of  a reflected 
binary  or  grey  code  corresponding  to  any  continuous  input 
intensity  value  between  0 and  8.  To  obtain  all  these  three 
, bits  in  parallel  we  can  place  an  array  of  3 periodically 


repeated  attenuating  strips,  which  can  attenuate  by  factors 
of  1,  i/2,  and  1/4  over  the  write  surface  of  the  LC  device. 

If  the  period  of  this  array  is  much  smaller  than  the 
inverse  of  the  maximum  spatial  frequency  of  the  input 
picture,  then  corresponding  to  each  input  intensity  value 
we  have  three  outputs.  These  outputs  are  the  three 
sinusoidal  curves  of  Fig.  6.3.  By  placing  an  array  of 
electronic  thresholding  devices  in  the  output  plane,  we  can 
sense  the  three  bits  at  the  same  time  all  in  parallel  as 
shown  in  Fig.  6.4.  Similar  ideas  have  been  used  for 
electro-optic  A/D  conversion  [6-2 , 8-3]  , but  none  have  been 
achieved  in  an  optically  controllable  device. 

6. 3 Analog- to-Digital  Conversion  Experiment 

In  this  experiment  we  want  to  obtain  the  bit  planes  of 
an  input  picture  in  real  time  as  explained  above.  The 
first  step  is  determination  of  the  response  curve  of  the 
liquid  crystal  device.  From  this  curve  we  obtain  suitable 
quantization  and  thresholding  levels  for  the  actual  A/D 
experiments. 

8.3.1  Device  Transfer  Function 

The  experimental  set  up  to  obtain  the  input  output 
characteristic  of  the  pure  birefringent  LC  device  is  shown 
in  Fig.  6.5.  The  input  light  source  is  a mercury  arc  lamp. 

A pair  of  fixed  and  rotatable  polarizers  in  the  input  light 
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Low  LCLV  Photodetector 

Spotiol  S Periodically  Array 

Frequency  Repeated  Electronic 

Input  Attenuating  Stripe  Threshold 


Figure  8.4  System  for  parallel  A/D  conversion 
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Figure  8.5  Experimental  set-up  for  real-time  parallel 
analog-to-digital  conversion 


beam  provide  a means  to  vary  the  input  light  intensity. 
The  Corning  short  wavelength  cut  filter  eliminates 
wavelengths  shorter  than  493  nm  to  make  sure  that  the  write 
beam  wavelength  is  within  the  sensitivity  range  of  the  CdS 
photoconductor.  The  read  light  source  is  another  mercury 
arc  lamp.  Due  to  the  dispersion  of  birefringence  in  the 
liquid  crystal,  it  should  have  a narrow  spectral  bandwidth. 
A cor  ion  interference  filter  with  peak  wavelength  434.7  nm 
and  bandwidth  of  18.4  nm  was  used  to  meet  this  requirement 
for  the  read  light.  With  no  picture  in  the  input  plane  the 
output  intensity  varies  in  a quasi-sinusoidal  fashion  with 
increasing  input  illumination  due  to  the  changing 
birefringence.  If  the  amount  of  birefringence  varied 
linearly  as  a function  of  the  write  beam  intensity,  one 
would  expect  a strictly  sinusoidal  variation  of  the  output 
intensity.  However,  a number  of  factors,  including  the 
optical  nature  of  the  liquid  crystal  and  the  photoconductor 
characteristic  properties  affect  the  output  characteristic 
and  produce  an  approximately  sinusoidal  output  whose 
frequency  varies  (monotonically)  with  input  intensity.  The 
response  curve  obtained  is  shown  in  Fig.  8.6. 

8.3.2  A/D  Experiment 

Although  the  theory  behind  the  A/D  conversion  assumes 
a strictly  periodic  response  characteristic,  it  is  in 
general  possible  to  produce  the  desired  bit  planes  using 

103 


00 


<u  g 
+-  o 
a 

E o 

2 g ® 

as  s 

<0  3 


Q .—  co 

•C  > 4> 

8 *o: 


\ « a 

\ o » 

\ -10 
\ 

\ 


•rH  a 

x>  m m 

a> 

<D  M XJ 
a)  -p 
M T3  -H 
x:  a>  » 

4J  M 
3 <U 
4)  tn  to 
x:  id  c 
4J  a)  o 
6 a 

M tn 

o a>  a> 
44  x:  >4 

4-> 

•o  Q) 

<u  tn  e 


a) 

<u  > a) 
u m x; 

•H  3 4-» 

> u 
a)  tn 

T3  T3  4-1 
•*4  C 
OH  # 

XI  o tn 
tn  tv 
ai  x 

sz  a)  a 
4-1  sz  <u 

H >4 


• ^ -H 
(I)  C M 4J 
> O 3 10 
4H0  3 
3 tn  C 
U WO  <V 
ai  a>  4-1 

0)  > 4-»  4-1 

tn  c 4->  a) 

coo 

O O TJ  V 

a ai 
in  Q u x 
<u  \x:  -h 
« < H 44 


A1ISN3Q  indino 


the  quasi-per iodic  response  curves  of  the  actual  devices. 
The  trade-off  that  must  be  made  is  that  one  must  resort  to 
a non-uniform  quantization.  The  quantization  levels  used  ' 
in  this  expriment  are  indicated  in  Fig.  6.6. 

There  are  no  attenuating  strips  on  the  LC  device  used 
in  this  experiment.  Instead,  the  bit  planes  were  generated 
serially.  Also,  the  output  was  recorded  on  hardclipping 
film  rather  than  a thresholding  detector  array.  A test 
target  was  generated  consisting  of  an  eight  gray  level  step 
tablet.  The  gray  levels  were  chosen  to  match  the 
quantization  levels  shown  in  Fig.  8.6.  This  test  object 
was  imaged  onto  the  LC  device  and  the  output  was 
photographically  hardclipped.  This  produced  the  least 
significant  bit  plane  of  a three-bit  A/D  conversion.  Next 
the  write  illumination  intensity  was  decreased,  effectively 
rescaling  the  device  response  curve  to  generate  the  next 
bit-plane.  The  last  bit-plane  (most  significant  bit)  was 
obtained  by  attenuating  the  write  intensity  again  and 
photographing  the  output.  The  input  and  three  bit  planes 
thus  generated  are  shown  in  Fig.  6.7. 

Although  the  output  contains  some  noise,  the 
experiment  illustrates  the  principle  of  real-time  parallel 
incoherent  optical  A/D  conversion.  It  was  found  later  that 
the  computer  generated  gray  scale  was  somewhat  noisy  due  to 
the  grain  of  the  high  contrast  film  used.  It  is  possible 


105 


that  future  experiments  with  cleaner  inputs  and  improved 
periodic  light  valves  should  produce  better  experimental 
results  and  more  bits  of  quantization. 

These  results  represent  the  first  real-time  parallel 
A/D  conversion  to  be  performed  on  two-dimensional  input. 
The  A/D  conversion  bit  rate  potential  can  be  estimated  from 
typical  parameters  of  devices  currently  available.  The 
important  parameters  are  device  resolution  which  is 
typically  40  line  pairs/mm,  device  size  which  is  on  the 
order  of  50mmx50mm,  and  speed  which  is  generally  designed 
for  30  frames/second.  Multiplying  all  the  parameters 

O 

together  implies  an  A/D  conversion  rate  of  1.2x10  points 
processed  per  second. 
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CHAPTER  9 


NONLINEAR  PROCESSING  USING  INTENSITY-TO-SPATIAL 
FREQUENCY  CONVERSION 

In  the  direct  method  of  nonlinear  optical  processing, 
described  in  the  previous  chapter,  the  halftone  screen  is 
eliminated  at  the  expense  of  losing  versatility.  In  this 
chapter  we  discuss  a method  of  real-time  nonlinear  optical 
processing  which  does  not  need  any  halftone  screen.  The 
method  operates  by  transforming  intensity  to  spatial 
frequency  before  performing  the  nonlinear  operation. 
Although  we  still  need  a coherent  optical  processor  for 
this  method,  its  advantages  of  eliminating  the  halftone 
screen  and  flexibility  makes  it  one  of  the  most  promising 
methods  for  real-time  nonlinear  optical  processing. 

In  the  following  sections  we  describe  a new  liquid 
crystal  device  capable  of  intensity-to-spatial  frequency 
transformation  and  describe  how  this  property  is  used  to 
obtain  an  overall  nonlinear  transformation  of  intensity. 
We  also  consider  the  processing  limitations  of  this  method 
and  present  experimental  results  showing  a real-time  level 
slice  nonlinearity. 
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9.1  Variable  Grating  Mode  (VGM)  Liquid  Crystal  Device 

In  chapter  two  it  was  mentioned  that  conduction 
effects  will  predominate  if  the  conductivity  of  the  liquid 
crystal  is  sufficient.  In  this  situation,  with  certain 
conditions  of  thickness,  voltage,  frequency,  etc.,  stable 
domains  are  formed  [9-1  to  9-3].  These  domains  are  like  a 
volume  phase  grating  structure.  Figure  9.1  shows  a picture 
of  such  domains  taken  through  a polarizing  microscope.  The 
spatial  period  of  these  domains  may  be  extremely  fine.  The 
period  of  the  grating  shown  in  Fig.  9.1  is  approximately  18 
ym.  The  interesting  property  of  these  domains  is  that 
their  period  is  a function  of  the  applied  voltage.  It  has 
been  observed  that  the  relationship  between  the  spatial 
frequency  of  these  domains  and  the  bias  voltage  is  linear 
over  a voltage  range  of  20  to  60  volts  and  spatial 
frequency  range  of  100  to  700  line  pairs  per  milimeter 
[9-3].  By  adding  a photoconductor  layer  to  the  liquid 
crystal  the  period  of  the  domains  should  be  a function  of 
the  light  intensity  falling  on  the  photoconductor.  The 
liquid  crystal  and  the  photoconductor  layers  can  be 
arranged  in  a form  similar  to  the  standard  LCLV  described 
in  Chapter  2 to  read  out  the  phase  grating  variation  in 
reflection.  Another  option  is  to  simply  sandwich  the 
photoconductor  and  the  liquid  crystal  layers  together  and 
use  different  wavelengths  for  the  write  and  read  beams. 
The  latter  configuration  has  been  employed  in  the  device 

109 

\ 


I 


Figure  9.1  VGM  domain  structure 
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used  in  our  experiment  as  shown  in  Fig.  9.2. 


9.2  Nonlinear  Processing  with  the  Liquid  Crystal  Variable 
Grating  Mode 

The  VGM  liquid  crystal  device  can  be  considered  to  be 
an  intensity-to-spatial  frequency  converter  capable  of 
operating  on  two-dimensional  images.  When  an  input  image 
is  imaged  onto  the  photoconductor  surface  of  this  device 
the  intensity  variations  of  the  input  image  change  the 
local  grating  frequency.  If  we  read  out  in  coherent  light 
and  Fourier  transform  the  read-out  image,  then  different 
spatial  frequency  components  of  the  encoded  image  would 
appear  at  different  locations  of  the  frequency  plane  as 
shown  in  Fig.  9.3a.  This  means  that  different  locations  in 
the  frequency  plane  correspond  to  different  intensity 
values  in  the  input  image.  Thus,  by  placing  appropriate 
spatial  filters  in  this  plane  we  should  be  able  to  obtain 
different  transformations  of  the  input  intensity  in  the 
output  plane  as  depicted  in  Fig.  9.3b.  This  figure 
describes  the  variable  grating  mode  nonlinear  processing 
graphically.  The  input  intensity  variation  is  converted  to 
spatial  frequency  variation  by  the  characteristic  function 
of  the  VGM  device  (upper  right-hand  quadrant).  This 
spectrum  is  modified  by  a filter  in  the  Fourier  plane 
(upper  left-hand  quadrant) . Finally,  the  intensity  is 
observed  in  the  output  plane  (lower  left-hand  quadrant). 
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Figure  9.3  Implementation  of  a nonlinearity 

utilizing  the  variable  grating  effect 
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Taken  together  these  yield  the  overall  nonlinearity  (lower 
right-hand  quadrant).  Design  of  a proper  spatial  filter 
for  a desired  transformation  is  a relatively  easy  task 
compared  to  the  design  of  a halftone  screen.  The  only 
parameter  available  in  the  design  of  the  spatial  filter  is 
the  location  in  the  frequency  plane.  For  example,  for  a 
level  slice  transformation  a simple  slit  that  passes 
certain  frequency  bands  will  be  adequate.  For  a better 
understanding  of  the  nonlinear  processing  with  VGM,  we  now 
formulate  this  process  in  mathematical  terms. 


Variable  Phase  Grating  Anal} 


Assume  that  we  have  a phase  grating  (periodic  array  of 
phase  shifting  elements)  extending  to  infinity  in  both  the 
x and  y directions.  Let  <j(x,y)  be  the  complex  amplitude 
transmittance  of  the  phase  grating.  Since  the  grating 
extends  to  infinity  we  can  assume  that  g(x,y)  is  a periodic 

function  of  x and  write  it  as  a Fourier  series  expansion 

00 

g(x,y)  = cexp(  j2imx/L)  (9,1) 

n=“°° 

where 

i rL 

-n  = L J 2<x'y>exP<-:>27rnx/L>dx  (9.2) 

and  L is  the  period  of  the  grating.  Because  we  assume  that 
the  grating  extends  to  infinity  in  the  y-direction,  £(x,y) 
is  not  an  explicit  function  of  y. 
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Now  we  assume  that  there  is  an  aperture  a(x,y)  in  the 
input  (x-y)  plane  where 


a (x,y)  = 


1,  for  (x,y)  inside  the 

aperture 

0,  otherwise 


(9.3) 


Then  the  transmittance  function  is 


+ CO 


t (x,  y)  = 2.(x,y) -a  (x  ,y)  = a (x,y)  • c exp  ( j2iTnx/L)  . (9.4) 

^ — n 


n=-c 


If  this  aperture  is  placed  in  the  input  plane  of  a coherent 


optial  processor,  then  in  the  Fourier  plane  we  get 

+oo 

— ( f X ^ f Y ) = A(fx,fy)* 


n=-°o 


(9.5) 


where 


T ( f x , fy)  = y(t  (x,y) } 

— ( f x ' f y>  = 


or 


+co 


T(fx,fy)  = 2 £nA(fx-g,fy) 

n=-°° 


(9.6) 


For  simplicity,  we  assume  that  the  aperture  is  square  with 
dimension  bxb.  Then 


a(x,y)  = rect (£) rect (£) 


(9.7) 


and 


A(fx,fy)  = b2sinc (bf x) sine (bfy) 


(9.8) 


Hence 


13.5 


-foo 

— ( f X ' f Y)  = ^ cnb2sinc[b(f  -£)]sinc(bf  ) 

n=-«>  Y (9.9) 

and  in  the  Fourier  plane  we  have  different  diffraction 

orders  each  of  which  is  a two-dimensional  sine  function. 


Note  that  if  b is  large  [9-4] 


and  we  have 


lim  b sine  bt  = 6 (t) 

b-*oo 

+oo 


(9.10) 


T(fx,fy) 


- ^£n5(fX-E)S(£Y> 

n=-oo 


-foo 


(9.11) 


n=-oo 


Assume  that  we  use  a filter  H in  the  Fourier  plane  that 
passes  the  first  order  only.  Then  the  amplitude  in  the 
output  plane  is 


V(x,y)  = y_1{c16(fx-i,fY)H(fx,fy)} 


- c^H  (jj,  0) 


The  output  intensity  is  then 


(9.12) 


VW  = 


(9.13) 


and  the  value  of  the  intensity  at  any  point  in  the  output 
plane  is  related  to  the  value  of  the  filter  function  at  the 
frequency  (i,0)  in  the  Fourier  plane.  Now  if  we  omit  the 
assumption  of  b->-<»  the  output  amplitude  is 
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(9.14) 


V(x0,y0)  = 3s  1{c1b2sinc(b(fx-i) ] •sinc(bfy) 

•H(fx,fy) } . 

This  assumes  that  the  orders  are  nonoverlapping,  a 
condition  which  occurs  if  b>>2L.  For  an  aperture  of 
general  shape 


V(x0,y0)  = y'1{c1A(fx-i,fy)H(fx,fy) } . (9.15) 

If  b is  the  average  dimension  of  the  apertur-;,  then  A will 
have  a width  on  the  order  of  2/b  as  in  the  case  of  the 
square  aperture.  The  output  intensity  in  this  case  is 

Io(x0»y0)  = |3t"1[c1A(fx-i,fy)H(fx,fy)]  |2.  (9.16) 

Now,  if  H ( f , f ) is  slowly  varying,  we  can  consider  it  to 

”*  A Y 

be  essentially  constant  over  the  width  of  A.  Thus  we  can 
write 

VVV  S |-'1[c1A(fx4,fy)H(i,0)]|2 

~ lcl^(fx_L'  fY^  ^ I 

1 2 (9‘17) 
= |H(i,0)c1a(x0,y0)exp(-j2TTx0/L)  | 

= c1|a(x(J,y0)  *H(i,0)  |2 

Going  back  to  the  definition  of  a(x,y)  in  Eq.  (6.3)  we  have 


To (x0'y0) 


cl|H(i,°) |2, 
0 


for  (xQ,y0)  within  the 
aperture 


otherwise. 


(9.18) 


I 
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Thus,  this  filtering  process  produces  an  aperture  in  ; 

the  output  plane  corresponding  to  an  aperture  a(x,y)  in  the 
input  plane  modulated  by  the  phase  grating  with  period  L.  , 

The  output  intensity  is  related  to  the  grating  period  i 

through  the  filter  function  given  by  Eq.  (9.18). 

The  above  analysis  can  be  used  to  estimate  the  , 

processing  limitations  of  the  VGM  device.  The  fundamental  > 

question  to  be  answered  is  how  large  the  smallest  picture 

element,  or  pixel,  must  be  with  respect  to  the  VGM  grating 
frequency. 

Suppose  that  the  usable  spatial  frequency  range  over 

which  the  VGM  device  can  operate  is  Av  . To  avoid 
_ o 

crosstalk  from  higher  orders  the  VGM  device  can  only  be 

operated  between  the  first  and  second  diffraction  orders  of 

the  lowest  frequency  used.  If  this  lowest  frequency  is 

denoted  by  \>  then  we  have 
o 

v = 2v  - v = v (9.19) 

o o o o 

as  shown  in  Fig.  9.4.  Now  if  the  number  of  intensity 

levels  that  we  wish  to  distinguish  is  N,  then  the  Fourier 

transform  of  any  one  aperture  represented  by  2 Av  in 

Fig.  9.4  must  be  contained  within  a region  of  width  v /N  in  ‘ 

o 

the  Fourier  plane.  So  we  must  have 

2Av  < vq/N  . (9.20) 
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Figure  9.4  Gray  level  resolution  with  VGM 


If  the  pixel  width  is  b then  the  width  of  its  Fourier 
tranform  is  2/b  which  means 


2Av  = | 
b 

Combining  Eq.  (9.20)  and  (9.21)  we  have 

B i V 


(9.21) 


(9.22) 


or 


bv  > 2N. 

o — 

This  requires  that  the  pixel  size  spans 
lowest  grating  frequency  if  N grey 
processed . 


(9.23) 

2N  periods  of  the 
levels  are  to  be 


As  an  example  if  =200  cycles/mm  and  N=50,  then 
each  pixel  must  have  a size  b=0.5  mm.  Thus  if  the  device 
has  a 50  mm  square  area  a 100x100  pixel  image  could  be 
processed . 


'j . 4 Experiment  with  the  VGM  Device 

In  this  experiment  we  show  the  ability  of  the  VGM 
device  to  generate  a level-slice  nonlinearity.  The 
experimental  setup  is  shown  in  Fig.  9.5.  The  input  picture 
is  illuminated  by  an  arc  lamp  source  and  imaged  onto  the 
photoconductor  surface  of  a VGM  device  which  has  an 
existing  phase  grating  structure  due  to  a bias  voltage. 
The  grating  period  is  locally  modulated  by  the  input 
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Figure  9.5  Experimental  set-up  for  real-time  parallel  level-slicing 


picture  intensity,  and  this  modulation  is  mapped  into  a 
position  along  a line  in  the  spatial  filter  plane.  A red 
filter  ensures  that  only  the  laser  beam  enters  the  coherent 
optical  processor.  Small  circular  annuli  of  varying  radii 
are  used  to  pass  certain  spatial  frequency  bands.  This  in 
effect  allows  certain  input  intensity  ranges  to  appear  in 
the  output.  Figure  9.6  shows  the  input  and  level  sliced 
output  pictures.  Figure  9.6a  shows  a positive  print  of  the 
original  image  as  photographed  on  the  imaging  screen.  A 
negative  of  the  original  was  used  in  the  experiments. 
Figure  9.6b  shows  a low  intensity  level  slice  corresponding 
to  a VGM  periodicity  of  120  lines/mm  with  approximately  3% 
bandwidth.  In  Fig.  9.6c  another  level,  corresponding  to 
153  lines/mm,  is  shown.  Figure  9.6d  at  236  lines/mm, 
illustrates  the  interference  from  second  harmonics.  Weak 
second  harmonics  of  the  low-intensity  image  slice 
corresponding  to  lib  lines/irm  can  appear  in  the  236 
lines/mm  level  slice.  In  Fig.  9.6e,  a broader  slice  of 
approximately  11%  bandwidth  was  taken  centered  about  the 
level  corresponding  to  140  lines/mm.  This  picture  may  be 
compared  with  the  previous  slices  and  particularly  with  the 
slice  shown  in  Fig.  9.6c.  Finally,  Fig.  9.6f  shows  a very 
high  intensity  slice  at  440  lines/mm  of  10%  bandwidth. 
Three  grey  levels  may  be  seen  simultaneously;  these 
correspond  to  the  superposition  of  three  wide  intensity 
slices. 
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Figure  9.6  Level  slicing  with  the  VGM  liquid  crystal  device 


CHAPTER  10 


CONCLUSIONS  AND  FUTURE  RESEARCH  TOPICS 

In  this  dissertation  we  have  presented  three  different 
techiques  for  performing  nonlinear  optical  processing  in 
real  time.  These  methods  specifically  use  the  following 
ideas  and  devices  to  obtain  nonlinear  transformations:  1) 
halftone  preprocessing  with  a real-time  image  transducer; 
2)  nonlinear  characteristics  of  a real-time  device;  and  3) 
capability  of  a real-time  device  to  produce 
intensity-to-spatial  frequency  conversion. 

To  implement  halftone  nonlinear  processing  in  real 
time  we  have  replaced  the  conventional  photographic 
hard-clipping  recording  medium  with  a real-time  image 
transducer.  The  nonbinary  characteristic  of  the  real-time 
device  degrades  the  output  of  the  process.  A detailed 
analysis  has  been  made  and  we  have  found  exact  methods  for 
predicting  such  degradations.  The  analysis  is  general  and 
works  for  any  recording  medium  characteristic  curve  shape 
and  any  halftone  screen  cell  shape.  It  has  been  shown  that 
for  some  monotonic  functions  it  is  possible  to  obtain  exact 
compensation  for  the  nonideal  characteristic  curve  of  the 
recording  medium.  The  compensation  problem  has  been  solved 
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in  general 

by  an 

approximate 

method 

which  obtains 

the 

halftone 

screen 

density 

profile 

by  minimizing 

the 

difference 

in  a 

mean- square 

sense 

between  desired 

and 

degraded  outputs.  The  advantages  of  the  halftone  method 
are  its  flexibility  to  produce  a variety  of  nonlinear 
transformations  and  that  it  is  a well-explored  technique. 
Its  limitations  are  the  soft  threshold  characteristic  of 
the  real-time  device  and  necessity  for  accurate  high 
space-bandwidth  product  halftone  screens.  A real-time 
logarithmic  transformation  via  this  method  has  been  shown. 

The  inherent  nonlinear  characteristics  of  a special 
real-time  birefringent  liquid  crystal  device  were  used  to 
obtain  direct  nonlinear  transformations.  This  method  is 
particularly  interesting  because  it  may  be  done  with 
incoherent  light  and  the  space-bandwidth  product  of  the 
image  to  be  processed  can  be  as  large  as  the 
space-bandwidth  product  of  the  device.  However,  it  is 
limited  because  the  type  of  nonlinearity  achievable  depends 
on  the  nonlinear  shape  of  the  characteristic  curve  of  the 
real-time  device  used.  A real-time  three-bit  parallel 
incoherent  analog-to-digital  conversion  nonlinearity  has 
been  obtained  by  this  method. 

The  property  of  the  liquid  crystals  to  generate,  under 
certain  conditions,  a volume  phase  grating  whose  period 
varies  with  the  applied  voltage  was  found  suitable  for 
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nonlinear  optical  processing.  It  was  shown  that  a 
real-time  device  with  such  property  is  capable  of 
intensity-to-spatial  frequency  conversion  which  can  be  used 
to  generate  nonlinear  transformations.  Like  the  halftone 
method  this  method  is  very  flexible  and  can  generate  many 
nonlinearities.  As  in  the  halftone  method,  the  carrier 
frequency  modulation  required  limits  the  space-bandwidth 
product  of  the  picture  to  be  processed.  A real-time 
level-slice  nonlinearity  was  obtained  by  this  method. 

Several  topics  seem  worthy  of  further  effort.  Our 
analysis  for  the  halftone  process  shows  that  although  we 
can  achieve  good  approximations  to  the  desired  nonlinear 
transformations  with  nonideal  devices,  a device  with  sharp 
threshold  characteristics  can  improve  the  results  quite 
significantly.  Research  for  real-time  image  transducers 
with  such  a characteristic  should  continue  so  that  the  full 
potential  of  the  halftone  method  can  be  achieved.  Finding 
new  methods  to  design  accurate,  high  space-bandwidth 
product  halftone  screens  is  still  a problem. 

The  capabilities  of  direct  nonlinear  processing  are 
not  yet  fully  understood.  New  applications  for  the 
existing  device  characteristics  remain  to  be  found  and  new 
devices  with  specified  characteristics  are  to  be  developed. 
One  obvious  extension  is  to  achieve  a parallel  A/D 
conversion  nonlinearity  with  more  than  three  bits.  In  this 
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1 


case  we  need  a device  which  has  four  or  more  intensity1 
peaks  with  very  regular  period  in  its  transfer 

I 

characteristic.  Other  parameters  like  wavelength  and  the 
bias  voltage  could  be  used  to  give  extra  dimensions  to  the 
whole  process.  For  example,  different  nonlinear 
characteristics  could  be  obtained  at  the  same  time  for 
different  wavelength.  The  immediate  use  of  such  a property 
is  in  A/D  conversion  where  we  then  can  obtain  different  bit 
planes  in  different  colors. 

In  the  area  of  the  VGM  device  the  theory  presented  in 
Chapter  9 should  be  extended  to  explictly  consider  entire 
images  as  a whole  rather  than  just  single  apertures.  This 
way  the  interaction  of  all  apertures  which  comprise  an 
image  will  be  taken  into  account.  Another  theoretical 
question  to  be  addressed  concerns  the  fact  that  the  VGM 
cdvices  do  not  produce  strictly  periodic  phase  gratings  as 
can  be  seen  in  Fig.  . This  results  in  a blurring  of  the 
diffraction  orders.  This  in  turn  affects  the  number  of 
gray  levels  that  can  be  processed  with  a given  device  and  a 
fixed  image  space-bandwidth  product.  The  three  dimensional 
nature  of  the  phase  grating  structure  observed  in  the  VGM 
device  is  another  factor  that  should  be  considered  in  a 
future  analysis  of  this  device. 
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APPENDIX  A 


To  simplify 

h (y) 

f (x) dx 

h(y-e) 

where  h ( • ) is  the  inverse  function  of  f(*),  let 
u = f(x)  =>  x = f_1(u)  = h(u) 

then 

x = h(y-e)  =>  u = f(x)  = f [h  (y-e ) ] = y-e 

and 

x = h(y)  ==>  u = f (x ) = f [h  (y)  ] = y 
also  du=f'(x)dx  but  since  x=h(f(x)) 

f(*>.hMf(*n  - 1 or  i*)  ■ inTrxrr = fptut 
==>  dx  = h' (u)du 

hence 

f My)  r y 

f(x)dx  = u*h' (u)du. 

•*  h(y-e)  * y-e 


J 

( A . 1 ) 


(A. 2) 


(A. 3) 


(A. 4) 


(A. 5) 


(A. 6) 
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t 


T 


Using  integration  by  parts  in  this  integral  we  have 


h (y) 


f (x) dx  = uh (u) 

J h(y-e) 


Y _ fy 

y-e  ■’y-e 


h (u)  du 


= yh (y) - (y-e ) h (y 


-■>-  f 

J \7  — 


h (u)  du 
y-e 


(A. 7) 
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